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Preface 


This  report  Is  the  result  of  my  attempt  to  design  a digital  control 
law  for  the  C-141  aircraft,  using  the  multivariable  time-domain  technique 
which  achieves  entire  elgenstructure  assignment.  Also,  an  Interactive 
computer  program  was  developed  that  aids  the  user  In  designing  the 
digital  control  law.  The  theoretical  background  that  I used  extensively 
came  from  the  work  of  Professor  John  D'Azzo  and  Professor  Brian  Porter  from 


the  University  of  Salford.  Two  previous  theses  by  Lt.  Stan  Larimer  and 
Capt.  James  Colgate  provided  assistance  In  the  development  of  the  Interactive 
computer  program. 

I wish  to  express  my  gradltude  to  Professor  D'Azzo  for  his  continuous 
guidance  In  the  development  of  this  thesis  and  his  painstaking  effort 
In  reading  the  thesis  for  accuracy  and  completeness.  I also  want  to 
thank  Major  Gary  Reid  and  Prof  C.  H.  Houpls  for  their  assistance  and 
suggestions. 


Thomas  A.  Kennedy 
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Abstract 

In  this  report,  a deadbeat  tracker  and  regulator  are  developed 
for  the  C-141  aircraft  for  three  different  flight  conditions  using 
the  method  of  entire  eigenstructure  assignment.  A heuristic  iterative 
method  is  presented  for  determining  an  "optimal"  sampling  time  for 
deadbeat  controllers  by  using  the  controllability  matrix. 

In  order  to  perform  the  above  designs  and  to  evaluate  them 
through  a conti nuous-time  simulation, a computer  aided  design  program  is 
developed.  This  program  is  fully  interactive  and  provides  complete 
error  protection  and  abort  protection.  The  program  allows  the  complete 
design  and  simulation  of  regulator,  tracker,  and  disturbance- rejection 
control  systems  with  full-order  observers.  In  addition,  the  program 
permits  complete  state  space  analysis  of  the  system  being  designed. 
Including  both  discrete  and  conti nuous-time  responses.  A user's 
manual  and  progranmer's  guide  are  provided  for  further  development  of 
the  program. 
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I.  Introduction 

The  increased  flexibilty  and  miniaturization  of  the  digital  computer 
has  made  it  a valuable  asset  in  the  design  of  modern  control  systems. 

In  order  to  make  full  use  of  the  digital  computer  as  a controller  for 
the  feedback  control  system, a discrete  control  law  methodology  must  be 
developed  and  proven  to  give  desirable  results.  A design  methodology 
that  meets  the  above  criteria  is  entire  eigenstructure  assignment.  This 
consists  of  the  assignment  of  the  closed-loop  eigenvalue  spectrum  and 
the  associated  set  of  eigenvectors  and  generalized  eigenvectors. 

Background 

Most  modal  control  approaches,  to  date,  deal  with  the  assignment  of 
eigenvalues  for  the  closed-loop  system  (Ref.  1).  The  ability  to  assign 
arbitrary  eigenvalues  to  the  closed-loop  system  allows  the  designer  to 
match  the  stability  characteristics  of  the  closed-loop  system  to  some 
desirable  model.  Therefore,  both  the  model  and  the  system  have  the  same 
steady-state  responses.  However,  there  Is  no  assurance  that  the 
transient  responses  of  the  closed-loop  plant  match  those  of  the  model. 

The  reason  for  this  Is  that  although  the  previous  modal  control  approach 
allows  arbitrary  assignment  of  an  eigenvalue  spectrum  to  match  some 
model.it  does  not  allow  any  flexibility  in  the  assignment  of  associated 
eigenvectors  of  the  closed-loop  system.  Thus,  the  closed-loop  system 
has  the  desired  stability  characteristics  but  does  not  necessarily  have 
the  same  transient  responses  or  performance  characteristics  of  the  model. 

Now,  with  the  entire  eigenstructure  assignment  methodology.  It  is 
possible  to  match  both  the  stability  and  the  performance  characteristics 
( of  the  closed-loop  plant  to  a desired  model.  Entire  eigenstructure 
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assignment  is  a time  domain,  multivariable  state  space  control  law 
methodology  that  assigns  both  the  desired  eigenvalue  spectrum  and 
associated  eigenvectors  to  the  closed-loop  system  (Ref.  2). 


Statement  of  Problem 

The  two  main  goals  of  this  investigation  are 

(1)  Develop  an  interactive,  user  oriented  computer  program  that 
aids  in  the  design  of  regulator,  tracker,  and  disturbance 
rejection  control  systems.  All  of  the  designs  are  performed 
in  the  time  domain  and  are  based  on  a multivariable  control 
law  methodology  utilizing  entire  eigenstructure  assignment. 

The  core  of  the  program  is  built  around  the  algorithms 
furnished  by  Professor  Porter  of  the  University  of  Salford, 

England  and  the  Air  Force  Flight  Dynamics  Lab  in  the 
computer  program  FORTRAC  (Ref.  3). 

(2)  Design  a regulator  and  tracker  for  the  longitudinal  axis  of 
the  C-141  for  three  different  flight  conditions,  using  the 
Interactive  computer  package  developed  in  goal  1.  All  of 
the  control  law  designs  are  based  on  the  C-141  having  both 
elevator  and  spoiler  control  surfaces.  The  method  of 
entire  eigenstructure  assignment  Is  used  and  all  closed-loop 
eigenvalues  are  assigned  to  the  origin  to  produce  deadbeat 

designs,  that  is,  the  output  response  reaches  steady-state  in  a finite 
number of  sampling  periods.  The  design  technique  presented  in  this 
thesis  also  produces  a ripple-free  response  In  the  steady-state 
for  constant  Inputs. 

Approach 

Structured  design  techniques  are  employed  stressing  modularity  and  top-down  design 
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in  the  formation  of  the  computer  program.  A requirements  definition  is 
first  developed,  stating  the  need  for  an  iterative  approach  to  the 
design  of  a control  law.  A detailed  explanation  of  the  structured 
design  approach  used  is  presented  in  Chapter  IV. 

In  the  design  of  the  control  laws,  the  C-141  Is  first  modeled  In 
the  conti nuous-time  domain  using  the  body  axes  stability  derivatives  found 
in  Ref.  4.  The  plant  is  then  discretized  using  a sampling  time  that 
is  determined  using  the  heuristic  approach  formulated  In  Chapter  II. 

Next,  for  the  tracker  design,  the  plant  is  augmented  with  discrete-time 
Integrators  to  produce  a zero  steady-state  error  in  the  desired  output. 

By  transforming  the  system  matrix  into  the  Brunovsky  canonical  form, 
the  control  Indices  of  the  system  are  determined  (Ref.  5).  This  Is 
done  so  that  the  minimum  possible  index  of  nilpotency  of  the  closed-loop 
plant  matrix  can  be  determined.  A deadbeat  control  law  is  then  synthesized 
by  assigning  all  eigenvalues  to  the  origin  In  order  to  drive  the  plant 
to  the  desired  states  in  p sampling  times,  where  p is  the  minimum 
Index  of  nilpotency  for  the  closed-loop  system  (Ref.  5).  Entire 
elgenstructure  assignment  Is  the  control  law  methodology  used  In  the 
control  law  synthesis.  Since  the  original  plant  is  a continuous  sytem, 
a conti nuous-time  simulation  Is  performed  using  the  sampled-data  control 
defined  above.  Thus,  the  simulation  depicts  the  conti nuous-time  response 
of  the  plant,  both  between  and  at  the  sampling  times.  Chapter  III  gives 
a detailed  explanation  of  the  above  design  methodology. 

Review  of  Literature 

Most  of  the  theory  presented  In  this  thesis  Is  from  works  by  Porter, 
Bradshaw,  and  D'Azzo.  Porter  and  Bradshaw's  work  on  nilpotency  properties 
in  the  design  of  deadbeat  controllers  Is  used  to  design  closed-loop 
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systems  that  have  a finite  settling  time  equal  to  the  size  of  the  minimum 
permissible  polynomial  (Ref.  5).  The  use  of  integral  plus  proportional 
feedback  In  the  design  of  multivariable  tracking  and  disturbance-rejection 
systems  by  Porter  and  Bradshaw  Is  drawn  upon  in  the  design  of  trackers 
that  have  zero  steady-state  error  in  the  desired  output  (Ref.  6). 

Also,  D'Azzo's  work  on  entire  eigenstructure  assignment  Is  used 
extensively  in  the  design  of  the  trackers  and  regulators  In  this  thesis 
(Ref.  2). 

Development  of  the  interactive  computer  program  is  based  on 
structured  design  approaches  recommended  by  Yordon  (Ref.  7). 

The  actual  Implementation  of  the  interactive  program  on  the  Wright- 
Patterson  CDC-6600/Cyber-74  was  aided  by  previous  work  by  Colgate  (Ref.  8) 
and  Larimer  (Ref.  9). 


II. 


In  order  to  design  a discrete  control  law,  the  continuous  plant  must 
be  discretized  using  an  appropriate  sampling  time,  T.  A theorem  that 
is  often  used  in  determining  an  appropriate  sampling  time  is  Shannon's 
Sampling  Theorem.  Shannon's  Theorem,  in  essence,  states  that  the 
sampling  frequency  &>s  must  satisfy  the  condition 


ois  > 2 uc  (1) 

where  o>$  is  the  sampling  frequency  and  ojc  is  the  highest  frequency 
contained  in  the  signal  to  be  sampled. 

However,  following  Shannon's  Sampling  Theorem  in  the  selection  of 
a sampling  time  to  discretize  a linear  multivariable  plant  does  not 
necessarily  guarantee  that  the  discretized  plant  Is  controllable  even  if 
the  original  continuous  plant  is  controllable.  Another  fact  that  must 
be  considered,  when  choosing  a sampling  time,  is  the  "degree  of  controllability" 
that  the  discretized  system  possesses.  In  other  words,  the  degree  of 
controllability  determines  the  magnitudes  of  the  feedback  gains  necessary 
to  drive  the  plant  to  the  desired  state.  Since  it  is  well  known  that 
high  feedback  gains  result  in  high  overshoots  in  the  transient  response, 
a system  that  has  a "low  degree  of  controllability"  Is  defined  In  this 
report  as  a system  that  requires  large  feedback  gains  to  drive  the 
plant  to  some  desired  state.  Conversely,  a system  that  requires  small 
feedback  gains  to  drive  the  plant  to  some  desired  state  Is  defined  here 
as  a system  with  a "high  degree  of  controllability." 


Determination  of  Degree  of  Controllability 

When  state-variable  feedback  Is  implemented  with  the  control  law 

u(kT)  * K x(kT)  (2) 

for  a discretized  plant  represented  by 

x(kT+T)  3 F x(kT)  + G u(kT),  (3) 

the  closed- loop  plant  is  represented  by 

x(kT+T)  3 (F  + GK)  x(kT)  (4) 

where  F is  the  discretized  nxn  plant  matrix,  G Is  the  discretized  nxm 
input  matrix,  and  K is  the  mxn  feedback  gain  matrix.  A necessary 
condition  for  the  above  system  to  be  controllable  Is  that  the  rank  of 
the  controllability  matrix  be  equal  to  n,  the  dimension  of  the  plant 
matrix,  where  the  controllability  matrix  is  defined  as 

Q 3 [G,FG....  Fn-1G]  (5) 

The  above  condition  determines  whether  a system  is  controllable  or 
uncontrollable.  However,  for  a controllable  system,  it  is  not  possible 
to  gain  any  insight  Into  the  degree  of  controllability  of  a system  from 
this  approach.  For  a controllable  system  the  controllability  matrix 
reduces  to  the  following  nxn  matrix  containing  n linearly  independent 
< vectors 
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where  9-|  >92 • • • > 9m  are  the  columns  of  G and  d^dg.-.-d^  are  the  control 
Indices  of  the  discretized  system  (Ref.  10:81).  The  difference  between 
the  two  representations  of  the  controllability  matrix,  Eq.  5 and  6,  Is 
that  in  the  representation  of  Eq.  5 the  controllability  matrix  is  square 
with  dimension  n only  when  G is  of  dimension  nxl.  But  In  the  representation 
of  Eq.  6,  the  controllability  matrix  is  always  square  with  order  n.  Thus, 

It  Is  possible  to  take  the  determinant  of  the  controllability  matrix  In 
the  form  of  Eq.  6.  Now  If  the  magnitude  of  the  determinant  of  the 
controllability  matrix,  as  defined  by  Eq.  6,  Is  small  the  system  Is 
nearly  uncontrollable.  This  condition  implies  that  one  or  more  of  the 
columns  of  the  controllability  matrix  are  almost  linearly  dependent. 

Also,  if  the  magnitude  of  the  determinant  is  small,  then  at  least  one  of 
the  eigenvalues  must  have  a small  magnitude. 

Therefore,  the  magnitudes  of  the  eigenvalues  and  size  of  the  determinant, 
for  the  matrix  of  Eq.  6,  possess  information  on  the  "degree"  of  Independence 
of  the  columns  of  the  controllability  matrix.  This  fact  Is  used  to 
determine  the  "degree  of  controllability"  of  a discretized  system. 

Thus,  if  a system  is  discretized  using  two  different  sampling  times,  and 
one  discretized  system  has  a larger  controllability  determinant  than  the 
other,  the  system  with  the  larger  determinant  Is  deemed  to  have  a 
higher  degree  of  controllability  than  the  system  with  the  smaller 
determinant. 
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Results 

The  C-141  model  for  medium  altitude  cruise  is  discretized  for  eight 
different  sampling  times  and  augmented  with  discrete  Integrators  as 
discussed  in  Chapter  VI.  The  determinant  of  the  controllability  matrix 
and  the  average  absolute  value  of  the  control  gains  are  calculated 
for  each  sampling  time.  These  values  are  listed  in  Table  I.  All 
runs  are  performed  with  no  modification  to  the  null  spaces  as  discussed 
in  Chapter  VI.  As  shown  in  Table  I the  results  suggest  that  using  the 
"degree  of  controllability"  to  determine  an  optimum  sampling  time  is 
more  of  a qualitative  method  than  an  absolute  method.  However,  this 
assumption  is  based  on  the  fact  that  the  smaller  the  feedback  gains  are, 
in  general,  the  smaller  the  transient  overshoots  are  for  a system. 

Entire  eigenstructure  assignment  allows  manipulation  of  the  null 
spaces,  and  the  changing  of  the  null  spaces  has  a direct  effect  on  the 
size  of  the  feedback  gains  in  the  algorithm  used  to  compute  the  feedback 
gains. 

It  is  a matter  of  experience  to  choose  a proper  sampling  time. 

In  this  case  it  appears  that  sampling  times  between  1.50  seconds  and  2.0 
seconds  should  lead  to  good  results.  But  it  is  not  possible  to  change 
the  null  spaces  of  the  systems  discretized  with  sampling  times  between 
0.1  seconds  and  1.25  seconds  to  get  feedback  gains  comparable  to  those 
for  the  other  sampling  times  of  1.5  to  2.0  seconds.  Thus  it  can  be 
concluded  that  the  determinant  of  the  controllability  matrix  Is  directly 

II  II 

related  to  the  degree  of  controllability  of  a system.  It  can  also  be  used 
as  a qualitative  method  in  determining  an  "optimal"  sampling  time. 
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Plant  Constraints 

Two  aspects  that  must  be  taken  into  consideration  in  choosing  an 
appropriate  sampling  time  for  a deadbeat  controller  are  the  size  of  the 
minimal  polynomial  of  the  closed-loop  system  and  the  physical  constraints 
of  the  plant.  The  size  of  the  minimal  polynomial  plays  a large  factor, 
since  the  index  of  nilpotency  is  equal  to  the  size  of  the  largest  Jordan 
block  which  can  be  selected  by  proper  assignment  of  the  eigenvector/generalized 
eigenvector  chains  to  the  closed-loop  eigenvalue  spectrum  (Ref.  2).  In  this 
case  the  transient  response  given  by 

x(kT)  = ( F+GK) Lx (0)  (7) 

becomes  zero  in  L sampling  times,  where  L is  equal  to  the  size  of  the 
minimal  polynomial.  Since  the  entire  elgenstructure  assignment 
methodology  allows  the  designer  to  choose  the  degree  of  the  closed-loop 
minimum  polynomial  to  be  greater  than  or  equal  to  the  largest  controllability 
index,  d1 , there  is  flexibility  in  choosing  the  "fastness"  or  the  time 
required  to  reach  steady-state  by  the  deadbeat  controller. 

In  order  to  decide  on  the  proper  degree  of  "fastness"  for  a deadbeat 
controller,  the  physical  constraints  of  the  plant  should  be  considered. 

For  example,  the  trackers  designed  in  this  thesis  for  the  C-141  are  to 
track  an  altitude  rate  step  input.  If  a sampling  time  of  1.0  seconds  Is 
used  and,  since  the  minimal  polynomial  of  the  closed-loop  system  Is 
designed  to  be  equal  to  4,  the  aircraft  would  track  the  command  Input  In 
a deadbeat  manner  In  4 seconds.  However,  since  there  Is  a delay  of  one 
sampling  period  before  the  aircraft  control  surfaces  receive  any  control 
Input,  the  aircraft  actually  only  takes  3 seconds  to  track  the  commanded 
Input  In  a deadbeat  manner,  from  the  time  of  actual  control  surface  deflection. 
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It  Is  obvious  that  this  response  is  much  too  fast  for  a transport  aircraft 
such  as  the  C-141.  There  are  two  possible  ways  to  circumvent  this  problem. 

One  way  is  to  increase  the  size  of  the  minimal  polynomial  of  the 
closed-loop  system.  The  other  solution  is  to  use  a larger  sampling  time. 

Summary 

The  following  steps  are  a heuristic  guide  for  choosing  a sampling  time 
for  a deadbeat  digital  controller: 

1.  Determine  a sampling  time  with  a "high  degree  of  controllability". 

2.  Check  the  lowest  degree  of  the  minimal  permissible  polynomial  of 

the  discretized  system.  This  is  equal  to  the  largest  control  index,  dj 

3.  Consider  the  physical  constraints  of  the  plant  and  decide  whether 

the  deadbeat  controller  is  too  fast  or  too  slow.  If  the  sampling  time 
does  not  satisfy  the  above,  go  back  to  step  one  and  pick  a new  sampling 
time  accordingly. 
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III.  Design  Methodology  Using  Entire  Eigenstructure  Assignment 

The  design  procedure  for  a digital  controller  consists  of  four  main 
steps: 

1.  Form  the  discrete-time  model  of  the  plant  from  the  continuous- 
time model. 

2.  Augment  the  discretized  plant  with  discrete  integrators  In  the 
case  of  a tracker  or  disturbance  rejector  design. 

3.  Perform  the  control  law  synthesis  in  the  discrete-time  domain. 

4.  Simulate  the  continuous-time  closed-loop  system  and  evaluate  the 
performance  of  the  digital  controller. 

These  steps  are  used  in  this  thesis  to  achieve  desirable  controllers 
for  both  regulator  and  tracker  designs  for  the  C-141  aircraft. 

Sampled-Data  Transformation 

The  first  step  in  the  design  of  a digital  controller  is  to  discretize 
the  continuous  plant  depicted  in  Fig.  la  and  described  by 


x(t) 


t 


y(t) 


A x(t)  + B u(t)  + R d(t) 
C x(t) 


to  form  the  sampled-data  system  shown  in  Fig.  lb  and  represented  by 


(8) 


i 


t 

I* 

it 


- 


x(kT+T)  - F(T)  x(kT)  + G(T)  u(kT)  + DR(T)  d(kT)  (9) 

y(kT)  - C x(kT) 

where  T Is  the  sampling  time.  This  transformation  can  be  accomplished 
by  using  Eqs.  10,  11,  and  12. 

O 
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Figure  la  Continuous  - Time  System 


Figure  lb  Discretized  System 


u 


F(T)  ■ e 


e Bdt 


i 

* 

DR(T)  = 


eAt  R dt 


where  eA  is  the  state  transition  matrix.  Since  the  system  is  assumed 
linear  and  time  Invariant,  the  B and  R matrices  are  constant  and  can  be 


taken  out  of  the  integral.  Thus,  all  that  is  left  is  the  determination  of 
eAT  and  the  | eAt  dt.  Two  methods  are  employed  in  the  solution  of  the 

Jo 


above  problem  in  this  thesis.  $ 

The  first  method  consists  of  diagonalizing  the  plant  matrix  as  shown 

by 

A - X AX"1  (13) 

' 

where  X Is  the  modal  matrix  for  the  plant.  Using  the  above  simplification, 
the  equations  for  the  discrete  transformation  reduce  to 


F(T)  - X eATX  _1 


G(T)  • Xt  fT  eAt  dt]  X _1  B (15] 

Jo 

DR(T)  - X [ [T  eAt  dt]X  _1  R (16] 

Jo 

Since  the  modal  matrix  is  constant  it  can  be  taken  outside  the  Integral 

The  computation  of  eAT  and  I eAt  dt  Is  a scaler  problem  where 

Jo 


e dt  is  a scaler  problem  where 


A.j(t=l,2,...,n)  are  the  distinct  eigenvalues  of  A which  appear  on  the 
diagonal  of  the  A matrix.  This  method  is  not  general  since  a system 
with  multiple  eigenvalues  produces  a singular  modal  matrix.  The 
above  method  is  used  to  discretize  all  continuous  plants  for  the  control 
law  synthesis.  However, a general  method  of  discretizing  the  continuous 
plant  containing  multiple  eigenvalues  employs  the  infinite  series  form 


.AT 


T + AT  + 


A2T2 

TT 


AnTn 

~rf! 


(17) 


AT 

where  I is  the  identity  matrix.  The  Integral  of  e is  then  calculated 
using  the  following  infinite  series 


IT  + 


AT‘ 

T T 


ft"-1!" 

nT 


(18) 


Both  series  converge  rapidly  for  a small  T and  are  easy  methods  for  digital 
implementation.  Due  to  time  limitations,  this  method  was  used  only  for 
the  continuous-time  simulation. 


Augmentation 

The  use  of  proportional  plus  integral  state  feedback  In  multivariable 
control  systems  results  in  zero  steady-state  error  for  both  tracking 
(Ref.  11)  and  disturbance-rejection  systems  (Ref.  12)  if  the  augmentented 
system  remains  controllable.  Therefore,  both  of  these  systems  are  augmented 
with  discrete  time  Integrators  described  by 

Z(kT  + T)  * Z(kT)  + T [e'(kT)]  (19) 

where  e'(kT)  is  the  error  vector  that  is  to  be  driven  to  zero  in  the 
steady -state.  This  error  is  represented  by  Eq.  20  for  the  disturbance 
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rejection  system. 


( 


e'(kT)  = M x(kT) 


(20) 


where  M is  a pxn  constant  matrix  with  p less  than  or  equal  to  the  number 
of  controls.  This  restriction  on  p is  made  so  that  the  augmented  system  is 
controllable  (Ref.  6).  Thus,  the  disturbance  rejection  system  is  described 
by 

— 

x(kT+T) 

z(kT+T) 

I 


F , 0 


TM  , I 


x(kT) 

z(kT) 


u(kT)  + 


DR 


d(kT)  (21) 


y(kT)  = [C,  0] 


x(kT) 

z(kT) 


It  is  obvious  that  any  combination  of  states  represented  by  Mx(kT)  is 
driven  to  zero  in  the  steady-state. 

The  error  for  the  tracking  problem  is 


e'(kT)  = M x(kT)  - v(kT) 


(22) 


where  v(kT)  is  the  piecewise  constant  conmand  vector.  The  augmented 
system  for  the  tracker  is  represented  by 
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w. 


v(kT)  (23) 


”x(kT+T) 

F , o' 

’x(kT)" 

’g" 

0 

s 

+ 

u(kT)  + 

z(kT+T) 

mm  % 

TM  , I 

z(kT) 
» «• 

0 

m « 

E 

■ w 

y(kT)  - [C,  0] 


x(kT) 

z(kT) 


where  E is  a pxp  diagonal  matrix  with  e^  a -T  and  p is  the  number  of 
commands  to  be  tracked.  The  vector  v(kT)  represents  the  piecewise 
constant  command  Inputs.  In  order  for  the  augmented  system  to  be  controllable, 
the  number  of  command  Inputs  to  be  tracked  must  be  less  than  or  equal  to 
the  number  of  controls  (Ref.  6). 


Control  Law  Synthesis 

The  purpose  of  entire  elgenstructure  assignment  Is  to  assign, 
using  state  feedback,  the  closed-loop  eigenvalue  spectrum 


d(F+GK)  = {Xr  X2 Xp}  (24) 

and  the  associated  set  of  eigenvectors 

(X”j » X2»*  • • »Xn)  (25) 

Since  the  eigenvalues  and  the  eigenvectors  for  the  closed-loop  system  are 
related  by  (Ref.  2:23) 

■ 0 (1  - 1,  2,  ....  n)  (26) 


[F-X^.G] 


“1 


where 


= kxt  (i  = l,2,..n)  (27) 

and 


(i  « 1.  2.  ...  , n)  (28) 


is  a vector  that  lies  in  the  null  space  of  the  matrix 


S(X1)  = [F-  ,G]  Cl  ■ 1.  2,  ....  n)  (29) 

the  feedback  matrix  K is  given  by  (Ref.  2:17) 

K = • • • » dip]  f xr  X2.  • • • * Xp]  (30) 


Although  no  restrictions  are  placed  on  the  desired  eigenvalue  spectrum, 
the  ker  S(X^)  does  place  a constraint  on  the  eigenvector  that  is 
associated  with  the  assigned  eigenvalues.  The  reason  for  this  is  that  the 
ker  S(x^)  specifies  the  null  space  within  which  the  eigenvector  must  lie. 
Also,  the  eigenvectors,  which  are  the  columns  of  the  modal  matrix  X, 
must  be  independent  since  the  Inverse  of  X is  used  in  the  calculation  of 
the  feedback  gain  matrix  as  shown  In  Eq.  30. 

In  the  case  of  multiple  eigenvalue  assignment,  as  in  the  design 
of  a deadbeat  controller,  it  may  be  necessary  to  generate  chains  of 
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generalized  eigenvectors  satisfying  the  relationship 

F v (mji»j)l 

rc  ^ t n Xi  _ . (m44-l  J) 


[F-X1 I ,G] 


(j  ~ 1»  2,  ...  i i = 1,  2,  . . . i p) 

This  generates  k..  strings  of  vectors  associated  with  the  eigenvalue 
(0 

X.,  where  x is  the  Jtth  vector  in  the  jth  string  which  is  of 

1 i 

length  m^.  since  the  ker  S(X^)  occupies  a subspace  of  size  equal  to  m, 
the  number  of  controls,  and  therefore  only  m linearily  independent 
eigenvectors  may  be  generated  from  this  null  space  (Ref  2:23,25). 

Equation  31  allows  the  generation  of  eigenvector/generalized  eigenvector 
chains  from  the  original  null  space  given  by  Eq.  29,  so  that  the  size  of 
the  largest  control  index,  d^,  of  the  open-loop  system  may  be  retained 
for  the  closed-loop  system.  This  is  possible  if  (1)  the  eigenvalue 
spectrum  of  the  closed-loop  system  is  assigned  to  the  origin,  (2)  the 
combined  length  of  the  eigenvector/generalized  eigenvector  chains  is  set 
equal  to  the  number  of  system  states,  and  (3)  the  largest  chain  length 
is  set  equal  to  d^.  In  this  case  the  outputs  are  equal  to  the  inputs 
after  a finite  number  of  sampling  times,  that  is,  the  system  is  deadbeat. 
This,  allows  a deadbeat  response  to  take  place  in  q sampling  times,  where 
q is  the  order  of  the  minimal  polynomial  of  the  closed-loop  system  and  is 
equal  to  the  largest  control  index  of  the  pair  [F,6], 

Regulator 

The  control  law  to  be  synthesized  using  entire  eigenstructure 
assignment  for  the  regulator  system  is 


. 


u(kT)  = K x(kT) 


where  K is  the  feedback  gain  matrix  that  is  calculated  using  Eq.  30. 
The  closed-loop  regulator  system  is  then  represented  by 


x(kT+T)  = (F+GK)  x(kT) 


y(kT) = C x(kT) 
and  is  shown  in  Fig.  2. 


Distrubance  Rejector 


Using  the  design  methodology  described  above,  with  the  plant  now 
augmented  with  discrete  Integrators,  the  control  vector  u(kT),  for  the 
disturbance  rejection  system  of  Fig.  3 with  v(t)  = 0,  becomes 


u(kT)  = K^kT)  + ^ z(kT) 


This  control  law  results  in  the  closed-loop  system  represented  by 


xCkT+T)  | [F+G^.FKgl  |x(kT) 


z( kT+T)  TM  , I z(kT) 


d(kT) 


y(kT)  - [C,  0] 


x(kT) 

z(kT) 


This  system  possesses  both  proportional  and  Integral  control. 
The  closed-loop  system  is  stable  if  the  desired  eigenvalue  spectrum 
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x(t)*Ax(t)+Bu(t) 


SAMPLER 


chosen  Is  Inside  the  unit  circle.  Therefore, 

lim  (z(kT+T)  - z(kT))  = 0 (36) 

and  as  a result,  in  the  steady-state 

TM  x(kT)  = 0 (37) 

Thus,  if  M Is  equal  to  the  output  matrix  C,  all  the  outputs,  y(kT)  * 

C x(kT),  are  driven  to  zero  in  the  steady- state. 

Tracker 

The  control  law  for  the  tracker  is  identical  to  the  control  law  for  the 
disturbance  rejector  as  described  by 

u(kT)  = K]  x(kT)  + Kg  z (kT)  (38) 


and  the  closed-loop  system  formed  by  this  control  law  is 


xCkT+T) 

F+GKr  GKg 

'x(kT)" 

“o’ 

a 

+ 

z(kT+T) 

™ * XP. 

z(kT) 

E 

• • 

V(kT) 


y(kT)  • [C,  0] 


x(kT) 


L2(kT)j 


(39) 


where  v(kT)  is  the  discretized  coimand  input  vector  and  E « [e^  ] is 


a pxp  diagonal  matrix  w1the^»-T.  Also,  M Is  a pxn  matrix, 

where  p is  less  than  or  equal  to  m,  the  number  of  controls,  and  n is  the 


number  of  plant  states.  When  the  closed-loop  eigenvalues  are  within  the 
unit  circle,  it  is  obvious  that 

lim  [z(kT+T)  - z(kT)]  = 0 (40) 

t-*» 


Thus, 

lim  [TM  x(kT)  + E v(kT)]  « 0 (41) 

t-*» 


and  therefore  the  desired  combination  of  states 

M x (kT)  (42) 

representing  some  of  the  outputs  will  track  the  conmand  inputs  with 
zero  steady-state  error.  Figure  3 represents  the  tracking  system. 

Continuous  Simulation 

Since  the  plant  Is  continuous.  It  Is  necessary  to  perform  a continuous- 
time simulation  of  the  closed-loop  system  in  order  to  evaluate  the  effect- 
iveness of  the  digital  control  law.  This  is  necessary  because  the  complete 
closed-loop  system  Is  partly  discrete  and  partly  continuous.  Therefore  the 
system  Is  simulated  using  the  continuous  state  space  representation  with 
piecewise  constant  control  Inputs. 

In  order  to  determine  the  continuous-time  solution  to  the  state 
space  representation  of  the  linear,  time  invariant  system,  a discrete 
approximation  to  the  exact  solution  Is  used.  This  method  Is  used  to 
discretize  the  plant,  as  explained  In  the  preceding  section  on  sampled 
data  transformation,  with  a sampling  time  that  Is  less  than  or  equal  to 
TSAMP/12,  where  TSAMP  is  the  sampling  time  of  the  digital  computer  used 


In  the  control  law.  The  state  and  output  responses  are  obtained  from 
one  sampling  time  to  the  next  with  the  control  Input  u(kT)  held  constant. 
Then  the  control  Input  Is  updated  at  the  next  sampling  time  and  the  process 
Is  continued.  Therefore,  It  Is  possible  to  determine  the  system  responses 
between  the  sampling  times. 
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IV.  Development  of  Initial  Program  Guidelines 


CESA,  which  stands  for  Complete  Eigenstructure  Assignment,  is  an 
interactive,  user  oriented  computer  program  that  uses  the  multivariable 
control  law  methodology  of  entire  eigenstructure  assignment  in 
the  design  of  the  following  three  different  control  systems: 
regulator,  tracker,  and  disturbance  rejector.  This  chapter  deals  with 
the  requirements  definition,  program  design  emphasis,  and  the  system 
constraints  that  were  faced  in  the  early  development  stages  of  the 
computer  program.  Chapter  V describes  the  actual  program  structure  of  CESA. 
Appendix  A consists  of  a User's  Manual  for  CESA  and  Appendix  B contains 
a Programmer's  Manual  to  CESA. 

Requirements  Definition 

In  order  to  exploit  the  advantages  of  the  structured  design 
techniques  formulated  by  Yordon  in  Ref.  7,  a fundamental  definition  of 
what  the  computer  program  Is  expected  to  accomplish  Is  needed.  This 
fundamental  definition  is  called  a requirements  definition  (Ref.  7). 

The  requirements  definition  for  the  computer  program  is 

1.  The  computer  package  must  be  Interactive  and  user  oriented  and 
allow  for  an  iterative  design  process. 

2.  The  program  must  be  able  to  recover  from  user  input  errors. 

3.  A process  must  exist  that  allows  the  user  to  store  and  update 
all  input  data  for  some  future  time. 

4.  Entire  eigenstructure  assignment  Is  the  control  law  design 
methodology  used  In  the  design  of  regulators,  trackers,  and 
disturbance  rejectors. 

5.  There  must  be  a method  available  that  discretizes  the  continuous 
plant. 


6.  The  plant  must  be  augmented  with  integrators  for  the  tracking 


7. 


8. 


9. 


10. 


11. 


and  disturbance  rejection  systems. 

The  user  must  have  complete  freedom  in  forming  the  eigenvectors 
from  linear  combinations  of  the  basis  vectors  in  the  appropriate 
null  space. 

The  controllability  and  observability  of  the  open-loop  system 
must  be  checked. 

The  ability  to  calculate  the  control  indices  and  the  observer 
indices  of  the  open-loop  system  must  be  available. 

A continuous- time  simulation  must  be  Included  In  order  to 
evaluate  the  control  law. 

A plotting  package  must  be  Included  so  that  the  state  response 
for  the  continuous-time  simulation  can  be  plotted. 


Using  the  above  design  requirements,  a structured  chart  was  derived 
as  shown  In  Fig.  4.  The  structured  chart  of  Fig.  4 Illustrates  the  modules 
of  the  program  and  the  Interconnection  between  the  modules. 


Program  Design  Emphasis 

The  manner  in  which  a program  is  physically  divided  Into  modules  can 
significantly  affect  the  structural  complexity  of  the  resulting  program. 
Two  measures  that  can  be  used  to  judge  the  goodness  of  the  design  are 
coupling  and  cohesion.  These  two  concepts  form  the  basis  of  structured 
design  theory. 

Coupling,  the  measure  of  the  strength  of  Interconnections  between  one 
moduie  .and  another,  must  be  kept  at  a minimum  to  assure  an  acceptable 
level  of  Independence  between  modules.  Thus,  If  a program  has  a low 
level  of  coupling  between  modules,  debugging  or  modification  of  one 
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module  does  not  affect  another  module.  The  level  of  maintainability, 
the  extent  to  which  a system  can  be  easily  corrected  when  bugs  are  discovered 
in  a program,  and  the  modifiability,  the  ability  of  the  system  to  be 
changed  or  enhanced  to  meet  the  needs  of  a user  during  a program's  life 

1 1 

time,  of  a program  are  both  increased  as  the  level  of  coupling  is  reduced 
(Ref  7:76-94).  Some  factors  that  influence  the  level  of  coupling  are 
the  type  of  connections,  the  complexity  of  the  interfaces  and  the  type 
of  data  that  flows  between  modules.  Also,  the  use  of  global  data, 
conmons,  has  a dramatic  effect  on  increasing  the  level  of  coupling. 

In  order  to  keep  the  level  of  coupling  low,  the  use  of  global  data 
is  kept  at  a minimum.  Also,  when  the  use  of  global  storage  is  necessary, 
labeled  commons  are  used  instead  of  blank  commons,  since  labeled  commons 
do  not  raise  the  level  of  coupling  as  high  as  blank  commons. 

Cohesion,  the  degree  of  functional  relatedness  of  processing 
elements  within  a single  module,  also  has  a direct  effect  on  the  maintainability 
and  modifiability  of  a program  (Ref  7:95,  125).  A high  level  of  cohesion 
results  in  a high  level  Of  maintainability  and  modifiability.  The  level  of 
cohesion  is  kept  high  within  a module,  by  putting  only  functionally  related 
processes  within  a module.  However,  when  memory  could  be  saved  at  the 
expense  of  a loss  of  cohesion,  the  saving  of  memory  outweighs  the  loss 
of  cohesion. 

A top-down  design  is  employed,  a design  strategy  in  which  the  major 
functions  of  a program  are  identified  and  their  implementation  expressed 
in  terms  of  lower-level  functions,  thus  stressing  modularity  and  acceptable 
levels  of  both  coupling  and  cohesion.  However,  as  shown  above,  it  was 
sometimes  necessary  to  make  tradeoffs  in  order  to  save  memory. 
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V. . Program  Structure 
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In  order  to  fit  CESA  in  the  60,000g  core  of  memory,  it  is  necessary 
to  use  special  techniques  to  save  memory  space.  Two  of  the  main  techniques 
used  to  conserve  memory  are  overlaying  and  mass  storage.  Also,  since  the 
program  is  completely  interactive,  the  user  interface  is  an  important 
part  of  the  program  structure.  The  purpose  of  this  chapter  is  to  discuss 
how  and  why  overlaying  and  mass  storage  are  used  and  to  explain  the  functions 
of  the  user  interface. 

Qygrl&yl_ng 

Overlaying  is  a technique  that  divides  a large  program  into  smaller 
programs  that  fit  in  the  memory  space  available.  Whenever  one  of  these 
smaller  programs  is  needed,  it  is  loaded  Into  central  memory  by  the 
executive  routine.  When  the  executive  routine  is  finished  with  a smaller 
program  it  unloads  It  and  then  loads  the  next  program  needed.  The  executive 
routine  Is  called  the  main  overlay  and  the  other  programs  are  called  the 
primary  overlays. 

The  primary  overlays  are  stored  on  a mass  storage  device  by  the 
system,  so  they  take  up  central  core  only  when  they  are  loaded  by  the  main 
overlay.  Thus,  as  long  as  the  combined  memory  needed  by  the  main  overlay 
and  any  primary  overlay  Is  less  than  or  equal  to  60,000g  the  total 
program  fits  on  the  Interactive  system. 

Blank  and  labeled  commons  are  used  to  pass  data  between  overlays. 

However,  to  keep  coupling  at  a low  level,  as  mentioned  in  Chapter  IV,  labeled 
comnons  are  used  Instead  of  blank  commons  for  communication  between 
overlays  In  CESA. 

A more  detailed  description  of  the  overlays  used  In  CESA  is  given  In  the 
Programner's  Manual  located  In  Appendix  B. 
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Data  Base  .1 

The  manner  in  which  a program  handles  data  plays  a strong  part  in  its 
efficiency  and  therefore  is  an  important  part  of  the  program  structure. 

Three  different  types  of  data  storage  are  used  in  CESA:  local  storage, 
global  storage,  and  mass  storage. 

The  two  different  types  of  mass  storage  used  are  sequential -access 
files  and  random-access  files.  Each  of  the  above  mentioned  storage  devices 
have  strong  points  and  weak  points  as  pointed  out  below: 

1.  Local  storage  is  temporary  storage  that  a program  uses  to  store 
variables  that  are  needed  only  during  the  execution  of  a particular 
routine.  Once  the  routine  Is  finished,  all  variables  stored  In 
local  storage  are  lost. 

2..  Global  storage  is  storage  that  each  routine  has  access  to.  On  the 
CDC  6600,  labeled  commons  and  blank  conmons  perform  this  function. 

Two  disadvantages  of  global  storage  are: 

a.  Memory  space  Is  taken  up  at  all  times,  whether  the  stored 
variables  are  used  or  not  used. 

b.  Global  storage  creates  a high  degree  of  coupling  between 
routines  as  mentioned  in  Chapter  IV. 

3.  Sequential -access  files  are  mass  storage  devices  that  data  may  be 
written  to  or  from  in  a sequential  format.  An  advantage  of 
a sequential  file  is  that  it  can  hold  much  more  Information  than 
global  storage  while  taking  up  the  same  amount  of  core. 

In  fact,  the  only  core  needed  for  a sequential  file  Is  Input-output 
buffer  core  (202g-2000g)  Ref  (13:16-2). 

A disadvantage  of  the  sequential  file  Is  that,  to  get  Information 
stored  on  the  file,  the  whole  file  must  be  read. 


32 


4.  Random-access  files  are  mass  storage  devices  which  differ  from 
sequential  files  In  that  Information  Is  written  to  a random-access 
file  In  a random  manner.  The  main  disadvantage  of  the  random  file 
Is  the  I/O  time  required  to  read  or  write  from  a random  access  file. 
Local  storage  is  used  In  CESA  as  temporary  scratch  registers  In  different 
routines.  Although  global  storage  takes  up  memory  and  increases  the 
coupling  between  modules.  It  Is  used  since  reading  and  writing  to  global 
storage  takes  much  less  time  than  reading  and  writing  to  mass  storage. 
However,  since  labeled  commons  increase  coupling  at  a lower  degree  than 
blank  commons,  they  are  used  instead  of  blank  commons. 

Three  sequential -access  files  are  used  in  CESA  as  storage  devices. 
Character  strings  are  stored  In  one  sequential  file  as  part  of  the  input 
data  verification  process  In  CESA.  Another  sequential  file  serves  the 
purpose  of  a backup  storage  device  for  CESA's  labeled  commons.  The  use  of 
this  file  also  allows  the  user  to  store  all  Input  data  for  some  future  use. 
This  cuts  down  on  data  input  time.  The  third  sequential  file  used  by 
CESA  Is  the  answer  file.  At^the  user's  command,  all  of  the  program  output 
Is  written  to  this  file,  which  then  can  be  disposed  to  the  system  printer 
after  termination  of  the  program. 

Random-access  files  serve  two  main  purposes  In  CESA.  One  purpose  Is 
to  satisfy  the  need  for  storage  of  the  numerous  matrices  that  are 
used  in  the  entire  eigenstructure  assignment  algorithm.  Twenty-nine 
two-dimensional  matrices,  ranging  In  size  from  14x14  to  33x23,  are 
stored  In  one  random  access  file.  The  second  purpose  that  the  random- 
access  file  serves  is  to  store  the  results  of  the  simulation  that  must 
be  made  available  to  both  the  printing  routine  and  the  plotting  routine. 
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User  Interface 


In  the  requirements  definition  in  Chapter  IV,  it  is  stated  that  an 
interactive  user  interface  is  needed  to  perform  program  control,  data 
manipulation,  and  error  protection  and  recovery.  The  manner  in  which 
these  functions  are  implemented  is  described  in  the  following  sections. 

Program  Control  Interface.  The  program  control  Interface  for  CESA 
is  modeled  after  the  program  control  Interface  that  is  used  in  TOTAL, 
an  Interactive  control  design  package  Ref.  (9).  The  user's  interface 
is  in  the  form  of  option  numbers  and  commands  with  which  the  user  controls 
the  program.  Options  are  selected  by  the  user  so  that  certain  computational 
functions  of  the  program  are  performed.  An  input  data  verification 
routine  determines  that  the  input  is  an  option  number  and  then  passes  data 
to  the  main  executive  which  then  selects  the  proper  routine  to  accomplish 
the  computation  selected.  Comnands  allow  the  user  to  set  mode  control 
switches.  For  example,  if  the  command  ANSWER,  ON  is  entered,  all  program 
output  is  written  to  a sequential -access  file  called  ANSWER  until  the 
comnand  ANSWER,  OFF  is  entered.  In  the  case  of  a command  input,  the 
input  data  verification  routine  deciphers  the  alphanumeric  input  and  then 
sets  a mode  switch  accordingly.  Appendix  B,  Programmer's  Manual,  gives 
a detailed  explanation  of  options  and  commands. 

Data  Manipulation.  The  data  manipulation  interface  allows  the  user  direct 
access  to  the  global  storage  of  CESA  which  contains  the  continuous  and 
discrete  versions  of  the  system  plant.  Either  of  these  two  system  discriptions 
may  be  printed  out  or  changed  upon  request.  In  order  to  have  one  of  the 
plant  matrices  printed  out,  the  name  of  the  matrix  is  typed  in  under  the 
option  mode.  Any  value  stored  In  a plant  matrix  can  be  changed  by  equating 
that  storage  location  to  a numerical  value. 


The  user  also  has  access  to  29  matrices  that  are  stored  in  a random  - 
access  file.  To  determine  the  contents  of  these  matrices, the  user  types 
the  name  of  the  matrix  to  be  printed.  For  a listing  of  these  matrices 
and  more  detailed  information  on  the  data  manipulation  interface, 
see  Append fcies  A and  B. 

Error  Protection  and  Recovers 

Three  levels  of  error  protection  and  recovery  are  provided  by  CESA.  The 
first  level  verifies  that  all  input  data  is  of  the  proper  format  and,  if 
it  is  not,  requires  the  user  to  reenter  data  from  the  last  correct  piece 
of  data.  In  the  case  of  numeric  input,  the  second  level  checks  that  the 
numeric  input*  is  within  the  proper  range  and,  if  itfsnot,  the  user  is 
required  to  enter  a proper  value.  Since  it  is  impossible  to  preceive  every 
possible  problem  that  a program  can  encounter  due  to  a user  input  error 
or  a program  error.  It  is  Impossible  to  completely  debug  a large  program 
without  It  being  in  full  use  for  some  trial  time.  In  order  to  spare  the 
user  any  significant  problems  that  may  be  encountered  by  a system  abort,  a 
CDC  system  recovery  routine  Is  Included  In  the  CESA  computer  package. 

Thus,  if  a mode  error  is  encountered,  the  program  recovers  and  no  data  Is 
lost.  This  last  level  is  a catch  all  and  Is  not  Intended  to  be  the  main 
error  protection  and  recovery  routine.  Rather  it  is  a device  that  Is  used 
to  prevent  any  unnecessary  hardship  to  the  user  due  to  a program  flaw. 

More  information  on  the  three  levels  of  error  protection  and  recovery  can  be  found 
In  Appendix  B,  the  programmer's  guide. 
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VI.  Deadbeat  Control  Law  Synthesis  Using  Entire  Eigenstructure  Assignment 


In  this  Chapter,  two  deadbeat  controllers,  one  a regulator  and  the 
other  a tracker,  are  developed  for  the  C-141  aircraft  for  three  different 
flight  conditions-landing,  low  altitude  cruise,  and  medium  altitude  cruise. 

The  strength  of  the  multivariable  design  methodology,  entire  elgensturcture 
assignment.  In  the  design  of  a deadbeat  controller  for  a multi-input, 
multi -output  system  Is  exhibited  by  the  addition  of  spoilers  and  speed 
control  to  the  C-141  aircraft.  All  of  the  estimated  spoiler  stability 
derivatives  and  all  of  the  longitudinal  stability  derivatives  are  taken  from 
Ref.  4.  The  mathematical  modeling  of  the  C-141  aircraft  is  done  in  the 
body  axes  since  all  aircraft  sensor  signals  are  generated  In  body  axes 
coordinates.  A continuous-time  simulation  is  performed  on  each  of  the 
closed-loop  systems  to  determine  the  effectiveness  of  the  deadbeat  controllers. 

C-141  Model 

The  longitudinal  equations  of  motion  used  for  the  C-141  model  are 
V * ^ w + ^ w + Mu  u + Mq  0 + 6fi  + M6sp6sp  (43) 
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All  of  the  stability  derivatives  in  the  equations  are  with  respect  to  the 
body  axes.  Since  Ref.  4 supplied  only  the  spoiler  stability  derivatives 
with  respect  to  the  stability  axes*  the  following  equations  from  Ref.  14  are 
used  to  compute  the  spoiler  stability  axes  derivatives  with  respect  to  the 
body  axes 


(Z4SP>b  = Z«SP  C0S  ac  * « 


The  transfer  functions  for  the  elevator  and  spoiler  actuator-servo 


combinations  are  modeled  as 
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and  the  transfer  function  representing  engine  lag  is 
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In  order  to  express  the  dynamics  of  the  C-141  aircraft  in  state 
variable  form,  the  following  state  equations  are  used 


where 


The  above  system  Is  discretized  with  a sampling  time  of  2 seconds,  which 


is  selected  according  to  the  approach  mentioned  in  Chapter  II,  and  the 
discretized  system  is  described  by 

x(2k+2)  » F x(2k)  + G u(2k)  (56) 

All  of  the  coefficients  that  are  used  in  the  modeling  of  the  C-141  aircraft 
for  the  three  different  flight  conditions  are  given  In  Table  II  (Ref  4:178, 
197). 

tequlator  Design 

A deadbeat  regulator  is  designed  for  the  longitudinal  axes  of  the  C-141 
aircraft  for  three  different  flight  conditions  - landing,  low  altitude 
cruise,  and  medium  altitude  cruise.  In  order  for  the  system  to  exhibit 
a ripple-free  response  in  a finite  number  of  sampling  periods,  all  of  the 
eigenvalues  of  the  discretized  closed-loop  system  are  assigned  to  the 
origin.  The  sampling  time  used  in  all  cases  is  2 seconds. 

Landing 

A regulator  is  designed  for  the  landing  flight  condition.  Substituting 
the  values  of  the  stability  derivatives  into  the  state  equations,  Eqs.  54 
and  55,  the  continuous  plant  matrix  becomes 


and  the  difference  equation  for  a sampling  time  of  2 seconds  Is 
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All  of  the  eigenvalues  of  the  closed-loop  system  are  assigned  to 
the  origin  so  that  the  flight  control  system  exhibits  a deadbeat 
response.  The  control  Indices  of  the  matrix  pair  (F,G)  are  the  set 
of  Integers  {3,2,2).  In  order  for  the  closed-loop  system  to  achieve 
a ripple- free  response  in  three  sampling  periods,  three  eigenvector/ 
generalized  eigenvector  chains  of  lengths  three,  two,  and  two  are 
generated.  Thus,  using  the  algorithm  from  Ref.  2 
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The  vectors  selected  from  the  null  space  of  SCO)  are 
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The  first  set  of  generalized  vectors,  generated  by  using  the  method 
described  in  Chapter  III  from  Ref. 2,  are 


F + GK 


- .3185,  .0007,  .0001  ,-  .4762,  .0217, -.0111,  .00007 

51.12  ,-.0241 .-.0187  , 45.41  ,-3.621  ,2.821  , .0662 

-5.310  , .0484,  .2133  ,-13.13  , .3398, -.3054, -.7307 

.3298, -.0006, -.00002,  .4176,-  .0229,  .0137,  .0003 

.6049, -.0049,.. 0022  , 1.471  ,-  .0409, -.0117, -.0054 
1.150  ,-.0089,  .0051  , 3.606  ,-  .0711 ,-.0356, -.0130 
-1.972  , .0193,  .06114,-5.687  , .1232,-. 0581 ,-.2118 


All  of  the  eigenvalues  of  the  closed-loop  system  are  located  at  the 
origin  and  the  size  of  the  minimal  polynomial  Is  three  as  required  by  the 
assignment  of  the  associated  elgenvector/generallzed  eigenvector  chains 
from  Eq.  62  through  Eq.  68.  A continuous-time  simulation  of  the  closed-loop 
system  Is  performed  and  responses  are  plotted  In  Fig.  5 for  the  Initial 
condition 





which  represents  a 0.1  radian  (5.72  degrees)  initial  condition  for  the 
pitch  angle.  The  maximum  change  of  angle  of  attack  is  0.0233  radian 
(1.34  deg)  and  the  maximum  change  In  velocity  along  the  body  vertical 
axis  is  -2.031  ft/sec.  The  maximum  control  values  are  a 0.1471  radian 
(8.4  degrees)  elevator  deflection,  a 0.3541  radian  (20.7  degrees) 
spoiler  deflection  and  a -0.5687  RPM  change  for  the  engines.  As  shown 
in  Fig.  5,  the  closed-loop  system  achieves  a ripple-free  steady-state 
response  in  three  sampling  periods.  Thus,  the  C-141  aircraft  exhibits 
good  transient  responses  and  the  ability  to  return  to  an  equilibrium 
state  in  a finite  length  of  time. 

Low  Altitude  Cruise 

A regulator  is  designed  In  this  case  for  the  low  altitude  cruise 
condition.  Equation  C-2  gives  the  continuous  plant  matrix  formed  from 
the  derivatives  from  Table  II  and  Eqs.  54  and  55  and  the  difference 
equation  for  a sampling  time  of  2 seconds  takes  the  form  of  Eq.  56. 

The  discretized  system  matrices  are  given  in  Appendix  C in  Eq.  C-3.  The 
control  invariants  for  this  system  are  the  set  of  Integers  (3,2,2)  and 
the  set  of  choosen  vector/generalized  vector  chains  are  given  by 
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The  resulting  feedback  gain  matrix  Is 
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- . 2686 , . 952E-3 ,- . 249E-3 ,-  .9896,1.0  ,-.0527,  .276E-3 
-1.541  ,.701 E-2,-. 136E-2,-6.494  , .5740,  .3208,  .1408E-2 
-.4302  ,.301 E-2,  .1942  ,-4.248  , .1238,  .1959, -.2532 


(75) 


With  the  resulting  closed-loop  matrix  found  In  Appendix  C,  Eq.  C-8. 

The  continuous  simulation  uses  the  same  Initial  conditions  as  given  In 
Eq.  72.  The  results  as  depicted  In  Fig.  6 show  that  the  closed-loop  system 
reaches  a steady-state  ripple-free  response  In  three  sampling  periods  (6  seconds) 
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with  a maximum  elevator  deflection  of  -.0989  radian  £5.67  degrees) 
and  a maximum  spoiler  deflection  of  -.6492  radian  (-37.2  degrees).  The 
maximum  change  in  engine  RPM  is  -0.417  RPM.  Thus,  the  closed-loop 
system  exhibits  a good  transient  response  with  the  ability  to  bring 
the  aircraft  back  to  an  equilibrium  condition  with  acceptable  control 
surface  deflections  and  engine  speed  changes. 

Medium  Altitude  Cruise 

The  stability  derivatives  for  the  medium  altitude  cruise  condition 
from  Table  II  and  Eqs.  54  and  55  are  used  to  set  up  the  state  equations 
for  the  aircraft.  The  resulting  continuous  plant  matrix  is  given  in 
Appendix  C,  Eq.  C-4.  This  system  is  discretized  with  a sampling  time  of 
2 seconds,  resulting  in  the  discrete  state  equation,  Eq  56.  The 
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The  resulting  control  gain  matrix  Is 
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The  closed-loop  matrix  using  the  feedback  gain  matrix  of  Eq.  78  Is  given 
In  Appendix  C,  Eq.  Cr9.  Using  the  same  Initial  conditions  as  In  Eq.  72, 
the  continuous -time  simulation  was  obtained  with  the  results  shown  in  Fig.  7. 

Figure  7 clearly  shows  that  the  closed-loop  system  reaches  a deadbeat 

response  in  3 sampling  periods  (6  seconds).  The  maximum  elevator  deflection 

is  -0.1131  radian  (-6.48  degrees)  and  the  maximum  spoiler  deflection  is 

-.8908  radian  (-51.03  degrees).  The  maximum  change  In  engine  RPM  Is 

-0.2912  RPM.  The  controlled  aircraft  has  an  acceptable  transient  response  and  the 

desired  ability  to  return  to  an  equilibrium  condition  with  reasonable  control  surface 

deflections. 

Regulator  Design  Conclusions 

In  all  three  regulator  designs,  entire  elgenstructure  assignment 
proves  to  be  a powerful  tool  In  the  design  of  time-optimal  regulators 
for  the  longitudinal  axis  of  the  C-141  aircraft.  Using  a sampling  time 
of  2 seconds  and  an  Initial  pitch  angle  of  0.1  radian  (5.73  degrees), 


viators  of  all  three  designs  return  the  aircraft  to  an  equilibrium 
condition  within  3 sampling  periods  (6  seconds).  This  result  is  possible 
since  the  length  of  the  eigenvector/generalized  eigenvector  chains  are 
set  equal  to  the  control  indices = (3,2,2).  Table  III  summarizes 
some  of  the  salient  results  of  the  continuous-time  simulation  for  the 
three  regulator  designs. 

Tracker  Design 

In  this  section  three  time-optimal  tracking  systems  are  designed  for 
the  longitudinal  axis  of  the  C-141  aircraft.  One  tracker  Is  designed 
for  each  of  the  three  flight  conditions-landing,  low  altitude  cruise, 
and  medium  altitude  cruise.  Discrete-time  integrators  and  comparators 
are  added  to  the  discretized  system,  as  discussed  in  Chapter  III,  so 
that  the  aircraft  can  track  a command  input  with  zero  steady-state 
error.  The  aircraft  is  modeled  by  Eqs.  43  thru  48. and  Is  set  up  In 
state-space  form  according  to  Eqs  54  and  55  using  the  stability 
derivatives  from  Table  II. 

The  plant  in  this  section  is  augmented  In  order  to  provide  tracking 
of  an  altitude  rate  command,  where  altitude  rate  Is  described  by 

h = -w  + UQ  9 (79) 

It  Is  also  desired  to  keep  the  forward  velocity  and  the  spoiler  deflection 
at  the  equilibrium  values  when  the  aircraft  reaches  the  steady-state 
flight  condition.  Therefore,  In  all  three  tracker  designs  the  discretized 
plant  is  augmented  with  three  discrete-time  Integrators  and  comparators 
as  defined  by  the  equation 


TABLE  III 


Regulator  Characteristics  for  the  Longitudinal  Dynamics  of  the  C-141  Aircraft 


for  an  Initial  Condition  of  0.1  radian  for  Pitch  Angle 


PARAMETER 

LANDING 

LOW  ALTITUDE 
CRUISE 

MEDIUM  ALTITUDE 
CRUISE 

l5el|Mx(rad' 

0.1471 

0.0989 

0.1131 

lsspLx(rad) 

0.3541 

0.6492 

0.8889 

lSrpmlmaxtrP"t> 

0.5687 

0.4170 

0.2912 

|e|mM(r,d> 

0.0476 

0.0716 

0.0905 

l“Lx(rad> 

0.0233 

0.0534 

0.0777 

luLx(ft/sec) 

2.031 

1.317 

1.244 

Settling  Time  (Sec) 

6 

6 

j 

6 
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(81 


where  the  outputs  to  be  tracked  are 


y^kT) 

0,-1  ,o,uo, 0,0,0 

y2(kT) 

= Mx(kT)  = 

0,0,1,  0 ,0,0,0 

y3(kT) 

0,0,0,  0 ,0,1,0 

- 

■ 

x(kT) 


(81 


and  the  command  Inputs  are 


V'l  (kT) 

fi(kT)  command 

v2(H) 

= 

u(kT)  command 

v3 ( kT) 

6sp0cT)  command 

with 

h(kT)  = 12  radians/second 
u(kT)  = 0 

VkT)=  0 


(82: 


(83) 

(84) 

(85) 


The  augmentation  of  Eq.  80  to  the  discretized  plant  equation  results  In 
the  composite  discretized  Input  and  output  equations 


x(kT+T)  F , 0 x(kT)  G 0 

- + u(kT)  + v(kT)  , 

z(kT+T)  TM,  I z(kT)  0 E 


y(kT)  * [C.O] 


x(kT) 


z(kT) 


and  now  the  augmented  plant.  Input,  and  output  matrices  are 


TM  , I 


C - [C.O] 


As  shown  In  Chapter  III,  the  eigenvalues  of  the  composite  system  can 
be  assigned  arbltrally  if  and  only  if  the  unaugmented  discretized 
plant  is  controllable  and  the  number  of  discrete- time  Integrators,  p. 
Is  less  than  or  equal  to  the  number  of  controls  m. 


Landing 

The  continuous  plant  matrix  for  the  landing  condition  Is  given  by 
Eq.  57.  The  discretized  form  of  the  state  equation  for  a sampling 
Interval  of  2 seconds  for  the  landing  condition  Is  given  by  Eqs.  58,  59,  and  60. 
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Using  Eqs  79  thru  90,  the  augmented  system  matrices  are 


- .1224.-.937E-3,  .955E-3 

52.01  ,-.0616  ,-.1399 

-24.79  , .2471  , .8741 

•5935.-.277E-2,  .109E-2 

0 ,0  ,0 

0 ,0  ,0 

0 ,0  ,0 

0 #5.0  ,0 

0 ,0  ,2 

0 ,0  ,0 


,-  .0351,  .807E-2 ,-  .0148  .-.282E-2, 0,0,0 

, 7.796  ,-4.010  , 3.615  , .5091  ,0,0,0 
,-61.10  , 1.623  ,-1.097  ,-3.210  ,0,0,0 


.9759,-  .0414  , .965E-2.-.249E-2, 0,0,0 


, 0 

, 0 

0 

ft 

,399.7 
, 0 
, 0 


, .206E-8,  0 ,0  ,0,0,0 

0 , .336E-3.0  ,0.0.0 


, 0 , 0 

, 0 , 0 

. 0 , 0 

, 0 , 0 


, .0183  ,0,0,0 
,0  ,1,0,0 
,0  ,0,1 ,0 
,2  ,0,0,1 


.0386,- 


.498E-2 


(92) 


r 


0,  0,  0 

0,  0,  0 

0,  0,  0 

0,  0,  0 

0,  0,  0 

0,  0,  0 

0,  0,  0 

-2,  0,  0 

0,-2,  0 

0,  0,-2 


(94) 


Since  the  pair  (F,(a)  is  controllable  and  it  is  desired  that  the 
closed-loop  system  exhibit  a time-optimal  response,  all  the  eigenvalues 
of  the  closed-loop  system  are  assigned  to  the  origin  by  implementing  the 
control  law 

u(kT)  = K]  x(kT)  + Kg  z(kT)  (95) 

The  control  Indices  of  the  pair  (F,(a)  are  the  set  of  integers  (4,3,3}. 
Therefore,  by  selecting  a set  of  eigenvector/generalized  eigenvector  chains  with 
lengths  equal  to  the  control  indices  , it  is  possible  to  achieve  a ripple-free 
steady-state  response  for  all  states  in  4 sampling  periods  (S  seconds).  The 
null  space  of  S(0)  » [F,E],  using  the  algorithm  from  Ref.  2 is. 
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— rr-  ' * *1'  :■  " — : 

- • • • •• 


” ■■  -"r 


and  the  resulting  closed-loop  matrix  is 


[F+GK] 


- .5321  , 

.447E-2,  .338E-3,-  1.535 

, .0353 

,-  .1250,  .902E-3 

6.310  , 

- .110  , .1597  ,-2.218 

,-.7890 

,-31.34  , .2940 

-14.72 

, .2719  ,-.6923  ,-73.99 

, .9388 

,-  7.02  ,-.1184 

.993E-2, 

.430E-2 , .938E-3,- 

1.005  .-.250E-2,- 

.1658,  .145E-2 

1.098  ,- 

.0155  .-.952E-4, 

4.343  ,-.0724  , 

.2213.-.618E-2 

1.133  ,- 

.0304  , .215E-2, 

7.822  ,-.0706  ,- 

.3122,-.842E-2 

.207  , 

•129E-2,  .1456  , 

.7285,  .0110  , 

1.270  , .2758 

■2.0  ,0  ,399.7  ,0 

0 ,2.0  , 0 ,0 


, 0 ,0 

, 0 ,0 


. 0 ,0 


. 2 .0 


.433E-3,-  .058 

.1433  ,-17.44 
-.3576  ,-  3.545 
•789E-3,-  .0915 

-.151E-2,  .1186 

-.744E-3,-  .1421 

.0319  , .6894 

0 , 0 

1 . 0 

0 . 1 


(108) 


All  of  the  eigenvalues  of  the  closed-loop  system  are  at  the  origin  and,  since 
the  lengths  of  the  eigenvector/generalized  eigenvector  chains  are  equal 
to  the  control  indices  of  the  pair  [F,G},  the  order  of  the  minimal  polynomial 
for  the  closed-loop  system  is  equal  to  4.  The  state  responses  of  the  continuous 
simulation,  with  a command  input  vector  as  described  by  Eq.  82  thru 
Eq.  85.,  are  plotted  in  Fig.  8.  The  altitude  rate  command  and  the  altitude 
rate  response  are  plotted  in  Fig.  9.  It  is  clearly  seen,  from  the  plots, 
that  the  closed-loop  system  tracks  the  command  input  vector  in  4 sampling 
periods  (8  seconds).  The  plots  also  show  that  the  maximum  elevator  deflection 
is  -.0946  radian  (-5.42  degrees)  and  the  maximum  spoiler  deflection  is 
-.2648  radian  (-15.17  degrees).  In  the  steady-state  the  change  in  forward 
velocity  (u)  is  zero  and  the  spoiler  control  surface  returns  to  its 
equilibrium  position.  Therefore,  the  designed  control  system  for  the 
C-141  aircraft  tracks  the  command  vector  with  desirable  state  responses  and 
control  surface  deflections. 

Low  Altitude  Cruise 

Equation  C-2  gives  the  continuous  plant  matrix  and  Eq.  C-5  gives  the 
discretized  form  of  the  state  equations  for  a sampling  interval  of 
2 seconds  for  the  low  altitude  cruise  condition.  The  augmented  system 
matrices  are  formed  using  Eqs.  80  thru  91  and  are  given  by  Eq.  C-6.  In 
order  for  the  closed- loop  system  to  exhibit  a ripple-free  response, 
the  eigenvalue  spectrum  for  the  closed-loop  system  is  assigned  to  the 
origin  using  the  control  law  described  by  Eq.  95.  The  control  Invariants 
of  the  pair  [F,?Tj  are  the  set  of  Integers  (4,3,3).  Since  a time-optimal 
system  response  Is  desired,  the  lengths  of  the  eigenvector/generalized 
eigenvector  chains  are  set  equal  to  the  control  indices..  The  vector/ 
generalized  vector  chains  selected  from  the  null  space  of  S(0),  using  the 
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algorithms  from  Ref.  2,  are 


- .0453  ,-23.38 

,-  . 1 50E+5 , 

- .971E+7,  .0125  , 

9.337  , 

1.0  , 1.0 

,1.0 

1.0  ,1.0 

1.0 

- .5  ,0 

, o 

0 , - . 5 , 

0 

.275E-3,  .184E- 

•2,  . 184E-2, 

.184E-2,  .918E-3, 

.184E-2, 

- .1263  ,-86.04 

, . 558E+5 , 

- . 360E+8 , .059  , 

34.94 

- .0419  ,-41.44 

,-  .274E+5, 

- .177E+8,  .0489  , 

17.48 

- .3996  , 33.78 

,-  .209E+5, 

•135E+8.-.4505  ,- 

•11.34 

1 

0 

w 

1 

o 

,-1.0 

-i.o  ,r.  o 

1.0 

o 

• 

•k 

o 

, 1.0 

1.0  ,1.0 

1.0 

.0837  , 82.96 

, . 549E+5, 

•355E+8.-.0978  ,- 

•35.06 

0 ,-  .1263 

,-86.04 

- .558E+5.0 

.059  , 

.141 E-4 ,-  .028 

,-32.26 

- .215E+5.-.164E-4, 

.0431  . 

.746E-2,-  1.037 

,-354 .8  , 

- . 230E+6 , .841 E-2, 

-.2473  , 

6063. 

, .0192  , 8.764 

,5601 

1.0 

V 
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• 

0 

w 
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o 

, -1.0 
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,-  . 5 , 0 

, o 

•184E-2 

, .184E-2,-  .184E- 

■2,  .184E-2 

.225E+5 

, .0435  , 31.96 
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, .113E-2,  15.08 

, . 1 02E+5 

-8432. 

,-  .4069  ,-10.01 

,-7781 

1.0 

«• 

o 
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o 

, o 

1.0 

, 1.0  , 1.0 

, 1.0 

-.222E+5 

,-  .225E-2.-30.16 

, .204E+5 

34.94 

, 0 , .0435 

, 31.96 

13.77 

,-  .378E-6,-  . 393E 

-2,  11.66 

145.8 

, .759E-2,-  .2278 

,135.0 

84 
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and  the  resulting  feedback  gain  matrix  is 
K = [KrK2]  = 


L 0.1)  v (2.1)  v (3.1)  y l4"’)  V V O'2’  V (3.2)  V (1  *3)  V (2.3)  v (3,: 

I A0  ^AQ  *AQ  *AQ  *A0  *AQ  *AQ  *AQ 


. 0284 , . 1 63E-2 ,- . 309E-3 .8049.-.894E-2,-  .2810,  .412E-4.-.973E-3, 
-.5904, .0156  .-.124E-2.-9.308  , .2309  ,-1.359  , .473E-3.-.622E02, 
-.9142r .016  , .4366  , 4.370  , .2740  , .8829, -.3072  , .0102  , 


- . 1 46E-3 , — . 1 432 
-.388E-3.-.6063 
. 988E-1 , .5791 


(110) 


and  the  closed-loop  matrix  formed  with  this  feedback  matrix  is  given  in 

Eq.  C-9.  All  of  the  eigenvalues  are  at  the  origin  and  the  closed-loop 

system  has  a Jordan  Canonical  form  with  one  block  of  size  4 and  two  blocks 

of  size  3.  Therefore,  the  order  of  the  minimal  polynomial  of  the  system  is 

equal  to  4.  The  state  responses  of  the  continuous  simulation,  with  a command  input 

vector  given  by  Eqs  82  thru  85,  are  plotted  in  Fig.  10.  The  altitude 

rate  command  input  and  the  aircraft's  altitude  rate  response  are  plotted 

in  Fig.  11.  Effective  tracking  of  the  command  input  vector  is  achieved 

in  4 sampling  periods  (8  seconds)  as  demonstrated  by  the  plots.  Also, 

the  plots  show  that  the  maximum  elevator  deflection  is  .0234  radian 
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(1.34  degrees)  and  the  maximum  spoiler  deflection  Is  .1492  radian 
(8.55  degrees).  In  the  steady-state  the  change  In  forward  veloctly 
(u)  Is  zero  and  the  spoiler  control  surface  returns  to  Its  equilibrium 
position.  Thus,  the  controlled  C-141  aircraft  has  the  ability  to 
track  the  command  Input  vector  with  desirable  transient  responses  and 
control  surface  deflections. 

Medium  Altitude  Cruise 

The  continuous  plant  matrix  for  the  medium  altitude  cruise  condition 
Is  given  by  Eq  C-4  and  the  sampled-data  state  equation  for  the  medium 
altitude  cruise  condition  Is  given  by  Eq.  C-5.  Equation  C-7  gives  the 
augmented  sampled-data  equations  that  are  formed  using  Eqs.  80  thru 
91.  Since  It  Is  desired  that  the  closed-loop  system  have  a deadbeat 
response,  the  closed-loop  eigenvalues  are  all  assigned  to  the  origin. 

The  control  Indices  of  the  pair  [F,(Tj  are  the  set  of  Integers  (4,3,3). 
In  order  to  achieve  a time-optimal  closedrloop  system  response,  the 
lengths  of  the  eigenvector/generalized  eigenvector  chains  are  set  equal 
to  the  control  invariants.  The  vector/generalized  vector  chains 
generated  from  the  null  space  of  S(0),  using  the  algorithms  from  Ref.  2, 
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Thus,  the  feedback  gain  matrix  is 

k * [Kr  y = 


(!,2)f(0  (2.2)  (3,2) 

0 o o 


<2-3> 

0 *0 


X 


0.1)  „ (2,1)  (3,1)  „ (4,1)  w (1,2)  „ (2,2)  „ (3,2)  „ (3,2)  „ (1,3)  „ (2,3) 

1 
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The  entire  eigenstructure  assignment  methodology  and  the  addition 
of  integral  control  to  the  closed-loop  system  is  a powerful  method  for 
the  design  of  digital  flight  controllers  as  demonstrated  by  the  three 
time-optimal  altitude  rate  tracker  designs  for  the  C-141  aircraft. 

In  all  three  designs  a sampling  time  of  2 seconds  is  used  and, 
since  the  lengths  of  the  eigenvector/generalized  eigenvectors 
chains  are  set  equal  to  the  control  indices  {4,3,3}  of  the  pair  [^,5], 
the  aircraft  tracks  the  command  input  vector  with  a deadbeat, ripple- 
free  response  in  4 sampling  periods  (8  seconds).  Table  IV  depicts 
some  of  the  salient  values  of  the  system  responses  of  the  continuous- 
time simulation  of  all  three  designs. 


TABLE  IV 


Tracker  Characteristics  for  the  Longitudinal  Dynamics  of  the  C-141  Aircraft 

for  the  Command  Inputs 
fi  = 12  ft/s ec,  u = 0 ft/sec 


PARAMETER 


|SeLx(rad> 

l%lmax(rad) 

|5rpJmax<r'""> 
lei™,  ( (rad/sec) 


l“Lx  (rad> 

lul ( ft/sec) 


l^e'ss  <rad) 

l«splss  (rad) 

|SrpJ5s  ^ 

1 9 1 ss  (rad/sec) 

lelss(rad) 

l°lss  (rad) 

lu|$s(ft/sec) 

Ihl $s  (ft/sec) 
Settling  Time  (sec) 


= 0 radian 


LANDING 
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CRUISE 


.0234 

.1492 


.2773 

.0T57 


.0118 


.1059 
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Mfeiilllila 


12.01 


.0220 


CRUISE 


.0284 


.1981 


.4039 

.0163 


.0144 


.4404 


.0255 


.2786 


.0181 
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VII.  Conclusions  and  Recommendations 


Conclusions 

The  combined  use  of  the  sampling  time  selection  procedure,  presented 
in  Chapter  II,  and  the  method  of  entire  eigenstructure  assignment  provides 
a powerful,  time  domain  approach  for  the  design  of  deadbeat  controllers 
for  multi-input,  multi-output  systems.  The  distinct  relationship 
between  system  controllability  and  sampling  time  selection  is  an  important 
factor  in  the  selection  of  a suitable  sampling  time  for  the  digital 
controller.  Also,  as  mentioned  in  Chapter  II,  system  constraints  must 
be  considered  in  the  design  of  deadbeat  controllers.  In  the  case  of 
tracker  designs,  the  use  of  discrete  integrators  and  comparators  forces 
the  closed-loop  system  to  track  a constant  command  input  with  a deadbeat, 
ripple-free  response. 

CESA,  the  interactive  computer  package  for  the  design  of  digital 
controllers  using  the  entire  eigenstructure  assignment  methodology,  meets 
all  of  the  design  requirements  mentioned  in  Chapter  III.  It  allows  the 
user  to  completely  analyze  the  structure  of  the  state  space  representation 
of  the  plant  and  also  to  analyze  the  discrete  control  law  for  the  closed- 
loop  system.  Since  the  program  is  totally  interactive,  the  user  can 
use  an  iterative  design  approach  in  the  selection  of  an  "optimum" 
control  law.  This  attribute  is  considerably  important,  for  it  Is  usually 
necessary  to  repeat  the  synthesis  procedure  for  a number  of  selected 
eigenvector/generalized  eigenvector  chains  and  to  evaluate  the  resulting 
control  law  thru  a continuous- time  simulation  which  is  a part  of  the  program. 

Both  the  regulator  and  tracker  designs  for  the  longitudinal  dynamics 
of  the  C-141  aircraft  are  time-optimal  with  all  control  values  within  system 
limits.  All  of  the  tracker  designs  track  the  command  Input  vector  In  a 


number  of  sampling  periods  with  zero  steady-state  error.  The  designs  were 

evaluated  using  a continuous-time  simulation,  where  the  plant  is  continuous  with 
a digital  computer  as  the  controller.  Plots  of  all  of  the  design  responses 
are  included  in  the  thesis. 

Recommendations 

The  use  of  entire  eigenstructure  assignment  and  the  proper  selection 
of  a sampling  time  for  the  design  of  deadbeat  regulators  and  trackers  are 
investigated  in  this  thesis.  Also,  an  interactive  computer  program  is 
developed  to  aid  in  the  design  procedure.  However,  there  is  much  more 
room  for  further  work  in  all  of  the  above  areas.  Some  suggestions 
for  further  work  in  these  areas  are 

1.  The  development  of  specific  guidelines  in  the  selection  of  the 
eigenvectors  and  generalized  eigenvectors  from  the  null  space  of 
S(X.). 

2.  The  modification  of  CESA  by  adding  routines  that  determine  the 
multivariable  zeros  of  a system.  This  will  allow  further  Investigation 
into  the  affects  of  sampling  time  on  system  controllability,  since 

the  selection  of  a sampling  time  directly  affects  the  position  of 
the  zeros  of  the  discretized  system.  Also,  In  order  to  have  a 
minimum  phase  control  law  the  zeros  of  the  discretized  system  must 
be  inside  the  unit  circle. 

3.  The  addition  of  routines  to  CESA  that  allow  the  assignment  of 
complex  eigenvalues,  using  entire  eigenstructure  assignment. 

These  routines  are  presently  available  In  updated  versions  of 
FORTRAC  (Ref.  3). 


4.  The  addition  of  routines  to  CESA  that  allow  for  the  design  of 
linear  multivariable  discrete-time  systems  incorporating 
error-actuated  dynamic  controllers  (Ref. .17).  These  routines  are 
presently  being  added  to  FORTRAC  (Ref.  3). 

5.  The  investigation  of  the  use  of  the  method  of  entire  eigenstructure 
assignment  in  the  rejection  of  random  disturbances,  such  as  a 

gust  rejection  analysis  for  an  aircraft. 

6.  The  investigation  of  tracking  a C*  input  as  discussed  in 
Appendix  D. 
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Overview  of  CESA 

O 

CESA,  which  stands  for  Complete  Eigenstructure  Assignment  is  an 
interactive  computer  program  that  was  created  to  aid  the  user  in  designing 
a state-space  control  law  for  a multi-input,  multi-output  system  using  the 
method  of  entire  eigenstructure  assignment.  The  following  overview  describes 
some  of  the  highlights  of  the  program. 


*CESA  aids  in  the  design  of  three  different  control  systems: 

•A  regulator  is  designed  using  the  method  of  entire 
eigenstructure  assignment.  A full -order  observer  may 
also  be  designed  if  the  plant  states  are  inaccessable. 

•A  disturbance  rejector  is  designed  by  augmenting  the 
plant  with  integrators  and  comparators  and  then  determining 
a control  law  for  the  augmented  system  using  the  method  of 
entire  eigenstructure  assignment.  If  the  plant  states  are 
inaccessable,  a full-order  observer  may  be  developed. 

•A  tracker  is  designed  by  augmenting  the  plant  with 
integrators  and  comparators  and  then  determining  the 
control  law  for  the  augmented  system  using  the  method 
of  entire  eigenstructure  assignment.  A full-order 
observer  may  be  designed  if  the  plant  states  are  in- 
accessable. 

♦All  of  the  control  law  designs  may  be  evaluated  by  a continuous- 
time simulation  of  the  closed-loop  system.  Calcomp  plots  of  the  system 
state  and  output  response  are  available. 

♦A  complete  state-space  analysis  of  a desired  system  may  be 
performed. 

♦System  matrices  can  be  entered  in  either  continuous  or  discrete  form. 

The  continuous  system  may  be  discretized  using  any  desired  sampling  time. 
♦All  input  data  can  be  saved  on  a memory  file  for  future  use.  This 
prevents  the  laborous  task  of  inputting  the  same  data  whenever  the  same 
system  is  under  study. 
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♦The  program  performs  complete  error  detection  and  diagnostics,  thus 
preventing  the  loss  of  input  data  due  to  abnormal  termination. 

1 . Introduction  to  CESA 

CESA,  which  stands  for  Complete  Eigenstructure  Assignment,  is  an 
interactive  computer  program  package  that  uses  the  method  of  entire  eigenstructure 
assignments  the  design  of  control  laws  for  multi-input,  multi-output  systems. 

The  computer  package  contains  sixty  options  and  some  special  commands  that  give 
the  user  an  iterative  approach  method  in  the  design  of  control  laws  for 
regulators,  disturbance  rejectors,  and  trackers.  The  state-space  representation 
of  the  system  can  be  entered  in  either  its  continuous  or  discretized  form. 

Once  the  system  model  is  entered,  a complete  state-space  analysis  of  the 
system  model  can  be  performed.  A conti nuous-time  simulation  is  available 
to  evaluate  the  designed  control  laws,  and  Calcomp  plots  of  the  system  state 
and  output  responses  are  available.  The  following  should  give  the  user 
all  the  necessary  information  for  optimal  use  of  the  computer  program  package. 

1 . 1 CESA'S  Input  Modes 

CESA  has  two  input  modes  in  which  it  requests  the  user  to  input  information: 
OPTION  and  DATA.  Each  mode  has  its  own  restrictions  on  allowable  inputs 
and  its  own  method  of  requesting  information. 

OPTION  Mode.  The  OPTIONmode  is  the  main  command  mode  of  CESA  and  its 
prompt  message  is 
0PTI0N> 

Once  the  program  is  in  this  mode,  the  user  Is  allowed  to  select  any  option, 
perform  any  command,  modify  any  system  matrix,  list  any  system  matrix,  or 
list  any  design  matrix.  This  mode  is  the  main  input  mode  of  CESA  since 
It  allows  the  user  complete  access  to  the  program  and  data  files. 
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Multiple  commands  or  options  may  he  typed  on  one  line.  However  all 
input  items  must  be  separated  hy  a hlank  or  a comma. 

DATA  Mode.  CESA  enters  into  the  DATA  mode  when  information  is  needed 
to  perform  an  option.  A particular  prompt  message  is  printed  out  for  the 
user  depending  upon  the  information  needed.  There  are  two  types  of  DATA 
modes.  In  one  type  of  data  mode  the  program  requests  a yes  or  no  reply, 
and  one  of  the  two  must  by  entered.  The  other  type  of  DATA  mode  requires 
a numerical  input  and,  unlike  the  other  DATA  mode,  a user  abort  to  the 
OPTION  mode  may  be  made  from  this  data  mode  by  entering  a "$". 


1.2  CESA  Options 

CESA  contains  sixty  options  that  allow  the  user  versatility  in  the 
design  and  analysis  of  control  laws  for  regulators,  disturbance  rejectors, 
and  trackers.  The  Options  are  grouped  together  into  four  sets  as  follows: 
0-15:  Matrix  Input  Options 
16-29:  State-Space  Analysis  Options 
30-45:  Control  Law  Design  Options 
46-59:  Simulation  Options 

The  first  option  in  each  group  lists  the  options  of  that  group. 

There  are  also  two  auxilary  options  used  in  CESA.  Option  60  updates  all 
data  in  a backup  memory  file  and  option  100  gives  a brief  description  of 
CESA.  A complete  list  of  the  CESA  options  is  given  at  the  end  of  this 
manual . 

1.3  CESA  Matrices 

CESA  distinguishes  between  two  different  groups  of  matrices.  One 
group,  the  inputted  system  matrices,  may  be  modified  or  listed  anytime  during 
the  program  execution.  However,  the  second  group,  the  resulting  design 
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matrices,  may  not  be  modified  but  they  can  be  listed  at  any  time  during 
the  program  execution.  Also,  only  the  inputted  system  matrices 
are  saved  on  a memory  file  that  can  be  fed  into  CESA's  data  banks  at 
a latter  time.  The  following  is  a brief  description  of  the  matrices 
in  CESA. 


Inputted  System  Matrices 


AMAT 

6MAT 

CMAT 

CDMAT 

RMAT 

ASDMAT 

BSDMAT 

RSOMAT 

EMAT 


Continuous  plant  matrix 
Continuous  control  input  matrix 
Output  matrix 
Integrated  Output  matrix 
Continuous  disturbance  input  matrix 
Sampled-data  plant  matrix 
Sampled-data  control  input  matrix 
Sampled-data  disturbance  input  matrix 
Sampled-data  command  input  matrix 


Resulting  Design  Matrices 

AUGAMAT  — / 
AUGBMAT  — t 
AUGRMAT  — / 
AUGEMAT  — / 
AUGCMAT  — / 
QMAT  — ] 


QTRNMAT 

QCHMAT 

NULLMAT 


CONLAW 

QOBMAT 

QOBCHAN 

OBNULSP 

OBVMAT 

OBWMAT 


Augmented  plant  matrix 
Augmented  control  input  matrix 
Augmented  disturbance  input  matrix 
Augmented  command  input  matrix 
Augmented  output  matrix 

Initial  matrix  set  up  to  determine  the  subspace  of 
ker  S(X^) 

Matrix  containing  the  null  spaces  of  S(X^) 

Matrix  containing  the  changed  null  space  vectors 
Matrix  containing  just  the  null  space  vectors  of  S(X^) 

Matrix  containing  the  selected  eigenvectors  and 
generalized  eigenvectors  [x^] 

Matrix  containing  the  remaining  part  of  the 
selected  vectors  and  generalized  vectors  from 
the  null  space  of  S(X^). 

Feedback  gain  matrix:  WMAT  * VMAT”^ 

Initial  matrix  set  up  to  determine  the  subspace  of 
ker  S(X.j ) for  the  observer  design 

Matrix  containing  the  changed  null  spaces  of  S(X.) 
for  the  observer 

Null  space  of  S(X^)  for  the  observer  design. 

Matrix  containing  the  selected  eigenvectors  and 
generalized  eigenvectors  for  the  observer  design 
Matrix  containing  the  remaining  part  of  the 
selected  vectors  and  generalized  vectors  from 
the  null  space  of  SfXj)  for  the  observer  design 


OBTRAN 

-- 

Matrix  containing  the  null  spaces  of  S(Xj 
for  the  observer  design  1 

ADCMAT 

[A+LC]  where  A is  the  plant  matrix,  L is  the 
o observer  feedback  gain  matrix,  and  C is  the 
output  matrix 

CLMAT 

-- 

Closed-Loop  Matrix 

ACOMP 

— 

Composite  plant  matrix 

BCOMP 

— 

Composite  control  input  matrix 

ECOMP 

— 

Composite  command  input  matrix 

RCOMP 

— 

Composite  disturbance  input  matrix 

CNLCOMP 

• — 

Composite  control  law  matrix 

2.  Complete  Description  of  CESA  Options 

In  order  to  exploit  CESA's  assets  In  an  optimum  manner,  it  is  necessary 
to  have  an  understanding  of  the  60  options  available  in  CESA.  This 
manual  is  intended  to  provide  all  of  the  information  needed  to  accomplish 
the  above.  Each  <5f  the  following  options  described  may  be  selected 
by  simply  typing  in  the  option  number  while  the  program  is  In  the 
option  mode. 

2.1  Matrix  Input  Options 

Since  a control  law  design  can  not  be  started  without  knowledge  of 
the  plant  model,  it  Is  necessary  to  enter  the  state-space  representation 
of  the  system  for  which  the  control  law  Is  to  be  designed.  Therefore, 
CESA's  matrix  input  options  are  an  Integral  part  of  the  program. 

CESA  allows  the  user  to  input  either  the  continuous  state-space 
representaion  of  the  system 

x(t)  « AMAT  x(t)  + BMAT  u(t)  + RMAT  d(t)  (A-l) 

y(t)  » CMAT  x(t) 
y(t)  * CDMAT  x(t) 

where 

AMAT  — Continuous  plant  matrix 

BMAT  — Continuous  control  input  matrix 
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RMAT  ~ Continuous  disturbance  input  matrix 

CMAT  — System  output  matrix 

CDMAT  — Integrated  system  output  matrix 
or  the  discretized  version  of  the  state-space  representation 

x(kT+T)  * ASDMAT  x(kT)  + BSDMAT  u(kT)  + RSDMAT  d(kT) 

y(kT)  = CMAT  x(kT) 
y(kT)  = CDMAT  x(kT) 


where 


ASDMAT 

BASMAT- 

RSDMAT 

CMAT 

CDMAT 


Sampled-data  plant  matrix 
Sampled-data  control  input  matrix 
Sampled-data  disturbance  Input  matrix 
System  output  matrix 
Integrated  output  matrix 


Option  0:  List  options 

This  option  gives  a list  of  all  the  matrix  input  options  (0-15). 

Option  1:  Recover  all  Data  from  MEMORY  File 

During  the  execution  of  CESA,  the  user  may  update  the  MEMORY  file 
with  new  data  (Option  60),  and  at  the  termination  of  the  program  the  MEMORY 
file  Is  filled  with  data  representing  continuous  and  discrete  versions  of 
the  system  under  study.  Option  1 allows  the  user  to  recover  all  the 
data  stored  In  this  file.  Thus,  the  user  can  terminate  CESA,  save  the 
MEMORY  file,  and  at  a latter  date  restart  the  program,  using  option  1 to 
recover  all  data  from  the  MEMORY  file. 


Options  2-6  Continuous  System  Matrix  Input  Options 


These  options  allow  the  user  to  input  AMAT,  BMAT,  RMAT,  CHAT  and 
CDMAT.  All  of  these  input  options  have  identical  formats.  For 
example,  if  the  user  wishes  to  enter  the  command  Input  matrix, 

BMAT,  he  would  use  option  3. 


OPTION  >3 

INPUT  OF  ( BMAT)  MATRIX! 

ENTER  MATRIX  SIZE!  ROUSrCOLUMNS>  3»2 
ENTER  3 ELEMENTS  PER  COLUMN! 

COLUMN  1 > 0» 1*0 
COLUMN  2 > 0»0.1 


COL  > 1 

ROW 

1 0. 

2 1*000 

3 0* 


0. 

0* 

1.000 


And  If  he  wants  to  enter  the  integrated  output  matrix,  CDMAT,  he  would 
use  option  6. 


OPTION  >6 

INPUT  OF  < CDMAT  ) MATRIX!  * 

ENTER  MATRIX  SIZE!  ROWS r COLUhNS>  2»3 
ENTER  3 ELEMENTS  PER  ROW! 

ROW  1 > l.lrO 
ROW  2 > O.OpI 

COL  >1  2 3 

ROW 

1 1.000  1.000  0. 

2 0*  0.  1.000 


Note  that  in  the  first  example  there  are  fewer  columns  than  rows,  and  the 
user  is  requested  to  enter  BMAT  by  columns.  In  the  second  example 
the  number  of  columns  is  greater  than  the  number  of  rows  and  the  user  is 
requested  to  enter  CDMAT  by  rows. 

Option  7 Help  User  Set  Up  Continuous  State-Space  Model  of  System 
This  option  is  used  to  input  all  the  matrices  needed  for  the 
continuous  state-space  representation.  It  first  requests  the  user  to 
input  the  number  of  states,  number  of  control  Inputs,  the  number  of 
disturbances  and  the  number  of  outputs.  Then  the  program  determines  the 
size  of  the  matrices  and  requests  the  values  for  the  elements  of  each  matrix 
as  in  options  2 to  6.  Also,  the  user  is  asked  if  there  are  any 
outputs  that  are  to  be  integrated  for  either  tracking  or  disturbance  rejection. 
And  if  so,  the  number  of  integrated  outputs  are  requested  by  the 
computer.  Then  the  user  is  asked  to  input  the  elements  of  CDMAT  as  In 
options  2 to  6.  The  following  is  a typical  example  of  option  7. 


OPTION  >7.  - 

THE  STATE  EQUATIONS  TO  BE  ENTERED  HAVE  THE  FORMS 

XDOT<T>-<  AMAT  >X<T>  + < BMAT  >U<T>  + < RMAT  >D(T> 

¥4T>«<  CHAT  >X(TX 

AND  (CDMAT)  * X<T)  ARE  THE  OUTPUTS  TO  BE  PROTECTED 
FROM  DISTURBANCES  OR  TO  TRACK  COMMAND  INPUTS 

AND  X IS  A VECTOR  OF  N STATE  VARIABLES 
U IS  A VECTOR  OF  M INPUTS 
D IS  A VECTOR  OF  K DISTURBANCES 
Y 18  A VECTOR  OF  L OUTPUTS 

ENTER  NO.  OF  STATES » INPUTS  * OUTPUTS  * DISTURBANCES  > 3t2t2t2 

ENTER  AMAT  WITH  3 ROUS  AND  3 COLUMNS. 

ENTER  3 ELEMENTS  PER  ROUS 
ROU,l  > OtlrO 
R0U>2  > OrOrl 
ROM.  3 > t-6 


' 


ENTER  BHAT  WITH  3 ROWS  AND  2 COLUMNS 

ENTER  3 ELEMENTS  PER  COLUMN! 

COLUMN  1 > 0,1,0 
COLUMN  2 > OfOfI 


ENTER  CMAT  WITH  2 ROWS  AND  3 COLUMNS 

ENTER  3 ELEMENTS  PER  ROW! 

ROW  1 > l»i'0 
R0I&2 


18  THERE  A DISTURBANCE  MATRIX »RMATv 
ENTER  RMAT  WITH  3 ROWS  AND  2 COLUMNS 

ENTER  3 ELEMENTS  PER  COLUMN! 

COLUMN  1 > 1*0*0 
COLUMN  2 > 

COL  > 1 2 

ROW 

1 1*000  0* 

2 0.  1.000 


n 


ARE  THERE  ANY  OUTPUTS  TO  BE  INTEGRATED  EITHER  FOR 

ri?ro RJJ£CJLREJECTI0N  or  COHHAND  INPUT  TRACKING 
ENTER  YES  OR  NO  > YES 

NUNBER  OF  CONTROLS  FOR  8YSTEH  EQUALS  2 

hust  “ less 

ENTER  THE  NUNBER  OF  INTEGRATED  OUTPUTS  > 2 

E*T£R  CDHAT  WITH  2 ROUS  AND  3 COLUMNS. 


ENTER  3 ELEMENTS  PER  ROW: 
ROU  1 >1,1,0 
RON  2 > ~071F7r 


COL  > 1 

ROU 

1 1.000 

2 0. 


1.000 

0. 


0. 

1.000 


THE  8TATE  SPACE  REPRESENTATION  18  COMPLETE. 


Options  8-10  Sampled-Oata  System  Matrix  Input  Options 


Options  8 through  10  allow  the  user  to  Input  ASDMAT,  BSOMAT  and 


RSDMAT.  All  of  these  options  have  the  Identical  formats  as  options  2-6. 


Option  11  Help  User  Set  Up  Sampled-Data  State-Space  Represenatlon 


of System 


Option  11  Is  Identical  to  option  7 except  that  In  this  case  the 


sampled-data  system  representation  Is  entered.  The  following  Is  a 


typical  example  of  option  11. 

OPTION  >11 


THE  EQUATIONS  TO  BE  ENTERED  HAVE  THE  FORMS 

X(K+1  > - ( ASDMAT )X(K)  + (BSDMAT)U(K)  ♦ ( RSDHA.T  ) D < K ) 
Y(K)  ■ ( CHAT  )X(K) 

AND  (CDMAT)  « X(K>  ARE  THE  OUTPUTS  TO  BE  PROTECTED 
FROM  DISTURBANCES  OR  TO  TRACK  COMMAND  INPUTS 


AND  X IS  A VECTOR  OF  N STATE  VARIABLES 
U IS  A VECTOR  OF  M INPUTS 
D IS  A VECTOR  OF  K DISTURBANCES 
Y IS  A VECTOR  OF  L OUTPUTS 


ENTER  NO.  OF  STATES* INPUTS* OUTPUTS* DISTURBANCES  > 3. 2*2.2 


■ 


ENTER  ASDMAT  UXTH  3 ROWS  AND  3 COLUMNS 


ENTER  3 ELEMENTS  PER  ROW: 
ROWtl  > «*474..4530..073S, 
ROWt-2  >iT.3774f-.4&ll'.012^ 
ROW?*  >fe:0723i-*S73S»-.  1334 

u i/f 

COL  > 1 2 

ROW 

• 1 .7474  .4330 

2 -.3774  -.4411 

3 -.7230E-01  -.3733 


3 

•7330E-01 
• 1210E-01 
-.1334 


ENTER  B8DNAT  WITH  3 ROWS  AND  2 COLUMNS. 


ENTER  3 ELEMENTS  PER  COLUMN t 
COLUMN  1 >;.*3241. .4530.-1.041 
COLUMN  2 >. 0421.. 0735. .0121 

CbL  > 1 2 


ROW 

• 1 .3241 

2 .4330 

3 -1.041 


•4210E-01 

•7330E-01 

•1210E-01 


Enter 


ENTER 
ROW  1 


COL 


ROW 

1 

2 


CHAT  WITH  2R0W8  AND 


3 ELEMENTS  PER  ROW l 
>, 1.1.0- 
>ML. 

> 1 2 


1.000  1.000 

0.  0. 


3 COLUMNS. 


3 

0. 

1.000 


IS  THERE  A DISTURBANCE  MATRIX. RMAT.— YES  OR  NOT  > YES  _ 
ENlfc*  R8DMAT  WITH  3 ROWS  AND  2 COLUMNS. 

ENTER  3 ELEMENTS  PER  COLUMN: 

COLUMN  1 >i,9141. -.2524. -.4410, 

COliUHN  2 > -♦  3241 . ♦ 4530  . -1 . 0410 


l If 

COL  > 

1 

2 

ROW 

1 

.9141 

.3241 

2 

.2524 

.4330 

3 

• 4410 

-1*041 
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ARE  THERE  ANY  OUTPUTS  TO  BE  INTEGRATED  FOR  DISTURBANCE 
REJECTION  OR  TO  TRACK  COHHAND  INPUTS? 

ENTER  YES  OR  NO  YES.1 
NUMBER  OF  CONTROLS  EQUALS  2 

NO.  OF  INTEGRATED  OUTPUTS  MUST  BE  LESS  THAN  OR  EQUAL 
TO  THE  NUMBER  OF  CONTROLS*  ENTER  THE  NUMBER  OF 

INTEGRATED  OUTPUTS > 2 

ENTER  CDMAT  WITH  2 ROUS  AND  3 COLUMNS* 

ENTER  3 ELEMENTS  PER  ROW* 

ROW  1 > 1*1*0 
ROW  2 > 0.0.1 

V? 

COL  > 1 . 2 3 

ROW 

1 1.000  1.000  0* 

2 0.  0*  1.000 


THE  STATE-SPACE  REPRESENTATION  IS  COMPLETE. 


Options  12-15  Reserve  Options 

These  options  have  been  set  aside  in  case  future  modification  of 
CESA  requires  the  Input  of  additional  matrices. 

2.2  State-Space  Analysis  Options 

In  order  to  design  a control  law  for  a system  or  to  evaluate  a system's 
responses,  it  Is  necessary  to  analyse  the  state-space  representation  of  the 
system.  Options  17  through  29  perform  different  operations  on  the 
state-space  matrices  of  the  system. 

Option  16  List  Options 

Option  16  lists  all  of  the  state-space  analysis  options. 

Option  17  Compute  Eiqenstructure  of  AMAT 

This  option  evaluates  the  eigenvalues  of  AMAT  and  determines  the 
modal  matrix  of  AMAT  using  Elspack  routines  (Ref.  14). 


* 


Option  18  Check  Controllability  of  Continuous  System 


The  controllability  of  the  pair  [AMAT,  BMAT]  is  checked  by 
this  option. 


Option  19  Check  Observability  of  Continuous  System 


This  option  checks  the  observability  of  the  pair  [AMAT,  CMATJ. 


Option  20  Compute  Sampled-Data  System 


Option  20  computes  the  sampled-data  representation  of  the  system  under 
study  for  some  user  inputted  sampling  time,  T.  The  procedure  uses  the  fact 


ASDMAT  = exp  (AMAT*T)  = XeATX_1 
BSDMAT  = (X  [T  eAT  dx  X-1 )*(BMAT) 


RSDMAT 


= (X  [T  eAT  dx  X'V(RMAT) 


(A-3) 

(A-4) 


(A-5) 


where  T is  the  sampling  time,  X is  the  modal  matrix  for  AMAT,  and  A is  a 

diagonal  matrix  containing  the  eigenvalues  of  AMAT  along  its  diagonal. 

NOTE:  Since  the  inverse  modal  matrix  is  used  to  compute  the  state 

AT 

transition  matrix  (e  ),  the  sampled-data  representation  of  a matrix  with 
multiple  eigenvalues  cannot  be  computed  using  this  option. 

Option  21  Compute  Eiqenstructure  of  ASDMAT 


This  option  evaluates  the  eigenvalues  of  ASDMAT  and  determines  the  modal 
matrix  of  ASDMAT  using  Eispack  routines  (Ref.  14). 

Option  22  Check  Controllability  of  the  Sampled-Data  System 


Option  22  checks  the  controllability  of  the  pair  [ASDMAT,  BSDMAT]. 


, ' T\*,-  ,-*S, 

■ . ..  - •>  • 


[Option  23  Check  Observability  of  the  Sampled-Data  System 

This  option  determines  whether  the  pair  [ASDMAT,  CMAT]  is  observable. 

Option  24  Compute  the  Control  Indices  of  the  Sampled-Data  System 
Option  24  computes  the  control  indices  of  the  pair  [ASDMAT,  BMAT] 
by  transforming  the  system  into  the  Brunovsky  cannonical  form  (Ref.  5), 
and  the  counting  the  size  of  the  system  subblocks.  These  are  equal  to  the 
control  indices  of  pair  [ASDMAT,  BMAT].  Also,  using  the  control  indices, 
the  regular  vectors  of  the  controllability  matrix  are  used  to  form  an 
nxn  "reduced"  controllability  matrix,  where  n is  the  number  of  states 
in  the  system  and  a regular  vector  is  a linearly  independent  vector.  The 
determinant  and  the  eigenvalues  of  the  "reduced"  controllability  matrix 
are  evaluated  to  determine  the  "level  of  controllability"  of  the  system. 

Option  25  Compute  the  Observer  Indices  of  the  Sampled-Data  Systems 
This  option  transforms  the  pair  [ASDMAT',  CMAT']  into  the  Brunovsky 
cannonical  form  and  then  counts  the  sizes  of  the  transformed  system's 
subblocks  whlchare  equal  to  the  observer  indices  of  the  system. 

Option  26  Compute  the  Controllability  of  the  Augmented  System 

The  controllability  of  the  pair  [AUGAMAT,  AUGBMAT]  is  computed  by  this 

option. 

j 

I I 

Option  27  Compute  the  Observability  of  the  Augmented  System 

This  option  evaluates  the  observability  of  the  pair  [AUGMAT',  AUGCMAT']. 


Option  28  Compute  the  Control  Indices  of  the  Aqumented  System 

This  option  is  identical  to  option  24  except  AUGAMAT  replaces  ASDMAT  and 
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AUGBMAT  replaces  BSDMAT. 


Option  29  Compute  the  Observer  Indices  of  the  Augmented  System 

This  option  is  identical  to  option  25  except  AUGAMAT  replaces  ASDMAT  and 
AUGCMAT  replaces  CMAT. 


2.3  Control  Law  Calculation  Options 

CESA  allows  the  user  to  design  three  different  control  laws:  regulator; 
disturbance  rejector;  and  tracker.  In  cases  where  the  states  are  inaccessable, 
full-order  observers  may  be  designed.  The  method  of  entire  eigenstructure 
assignment  is  used  in  all  of  the  control  law  developments.  Therefore,  the 
user  can  assign  both  the  eigenvalue  spectrum  and  the  associated  eigenvector/ 
generalized  eigenvectors  to  the  closed-loop  system. 

Options  31  thru  36  aid  the  user  in  developing  continuous  control  laws 
for  the  above  systems.  Options  38  thru  43  allow  the  user  to 
design  discrete  control  laws  for  regulator,  disturbance  rejector  and 
tracker  systems. 


Option  30  List  Options 

This  option  lists  options  30  thru  45. 


Option  31  Design  Continuous  Regulator 

♦NOTE:  This  option  has  not  been  completed  at  this  writing. 


Option  32  Design  Continuous  Regulator  with  Full-OMer  Observer 
♦NOTE:  This  option  has  not  been  completed  at  this  writing. 


Option  33  Design  Continuous  Disturbance  Rejector 


•NOTE:  This  option  has  not  been  completed  at  this  writing. 


Option  34  Design  Continuous  Disturbance  Rejector  with  Full-Order  Observer 


♦NOTE:  This  option  has  not  been  completed  at  this  writing. 


(j 


Option  35 
♦NOTE; 


Design  Continuous  Tracker 

- H 

This  option  has  not  been  completed  at  this  writing. 


! 


Option  36  Design  Continuous  Tracker  with  Full-Order  Observer 
♦NOTE;  This  option  has  not  been  completed  at  this  writing. 

Option  37  Reserve  Option 

Option  38  Design  Sampled-Data  Regulator 

This  option  uses  the  entire  eigenstructure  assignment  methodology  in 
the  design  of  the  digital  control  law 

u(kT)  = K1  x(kT)  (A-6) 

where  T is  the  sampling  time,  and  K.  is  the  state  feedback  matrix. 

The  program  accomplishes  this  in  the  following  five  steps: 


O 


1.  Discretizes  the  continuous  system,  if  needed,  according  to  a user 
Inputted  sampling  time,  T. 

2.  Performs  a complete  state-space  analysis  on  the  system  upon  user 
request. 

3.  Allows  the  user  to  assign  distinct  real  eigenvalues  or  multiple 
real  eigenvalues  and  the  associated  eigenvector/generalized  eigenvectors 
to  the  closed- loop  system. 

4.  Determines  the  feedback  gain  matrix  K-|  according  to 

K1  " Cx-|»X2»***»Xn]  1 (A-7) 

* fl  X"1 

where 
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r •• • • T *>+  t yZ " 


e ker  [ASDMAT  - A^I,  BSDMAT]  (A-8) 

5.  Forms  the  closed-loop  matrix  CLMAT  where 

CLMAT  = [ASDMAT  + BSDMAT*^]  (A-9) 

Option  39  Design  Sampled-Data  Regulator  with  Full -Order  Observer 

This  option  allows  the  user  tc  design  a digital  controller  for  a regulator 
when  the  system  states  are  iinaccessable  and  the  system  is  observable. 

Thus,  for  a linear  multivariable  sampled-data  system  governed  by  the 
discrete-time  representation 

x(kT+T)  = ASDMAT  x(kT}  + BSDMAT  u(kT)  (A-10) 

and 

y(kT)  - CMAT  x(kT)  (A-ll) 

with  the  states  x(kT)  observable  in  the  output  y(kT),  it  is  possible  to 
design  a discrete-time  observer  which  satisfies  the  matrix  difference 
equation 

x(kT+T)  = ASDMAT  x(kT)  + BSDMAT  u(kT)  + L[Cx(kT)-y(kT)]  (A-12) 

/\ 

where  x(kT)  Is  the  observer  state  vector  (Ref.  5).  Subtracting  Eq.  A-12 
from  Eq.  A-10,  using  Eq.  A-ll,  and  designating  the  observer  error  states 
x(kT)  by 

x(kT)  =■  x(kT)  - x(kT)  (A-13) 

yields  the  observer  error  state  equation 
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x(kT+T)  = (ASDMAT  + L*CMAT)x(kT) 


(A-14) 


If  the  L matrix  is  chosen  so  that  the  eigenvalues  of  the  matrix  [ASDMAT+L*CMAT] 


are  assigned  within  the  unit  circle,  then 


Tim  xCkT)  = lim  (ASDMAT + L*CMAT)K  = 0 


(A-15) 


for  any  initial  condition. 


Therefore,  as  k-*=°,  the  observer  states  are  equal  to  the  sampled-data  plant 


states  and  the  observer  states  x(kT)  may  be  used  instead  of  the  inaccessable 
plant  states  x(kT)  in  the  control  law 


u(kT)  = k x(kT) 


(A-16) 


The  composite  closed-loop  system  with  the  discrete-time  observer  and  the 


control  law  of  Eq.  A-16  is  represented  by 


x(kT+T) 


ASDMAT+BSDMAT*K1 ,-BS0MAT*K 


x(kT+T) 


,ASDMAT+L*CMAT 


*1 

1 

J 


x(kT) 


x(kT) 


(A-17) 


The  program  allows  the  user  to  determi ne both  the  matrix  and  the  L 
matrix  according  to  the  procedure  of  option  38.  Then  it  forms  the  composite 


system  given  by  Eq.  A-17  where 


CIMAT 


ASDMAT  + BSDMAT*K1 ,-BSDMAT*^ 


,ASDMAT+L*CMAT 


(A-18) 


■ 


Option  40  Design  Sampled-Data  Disturbance  Rejector 

This  option  allows  the  user  to  design  a digital  control  law  of  the 


u(kT)  = K1  x(kT)  + K2  z(kT) 


(A-19) 


for  the  augmented  discrete-time  system 


x(kT+T) 

z(kT+T) 


ASDMAT,  0 
TM,  I„ 


x(kT) 

z(kT) 


BSDMAT  RSDMAT 

u(kT)  + d ( kT ) 

0 0 


(A-20) 


where  z(kT)  are  the  discrete-time  integrators  that  are  augmented  to  the 
system  and  the  p chosen  outputs  are 


M x(kT) 


(A-20) 


where  p must  be  less  than  or  equal  to  m,  the  number  of  controls,  in  order 
to  maintain  controllability.  This  system  rejects  the  piecewise-constant 
disturbance  input  vector  d(kT)  (Ref.  11).  Once  the  augmented  system  is 
set  up,  the  program  procedes  as  in  option  38  to  form  the  closed- loop 


matrix 


CLMAT 


ASDMAT  + 3SDMAT*K1 , BSDMAT  * K? 


. I. 


(A-21) 


* * • \. 


Option  41  Design  Sampled-Data  Disturbance  Rejector  with  Full-Order 
Observer 

This  option  allows  the  user  to  design  a control  law  that  rejects 
plecewise-constant  disturbance  Inputs,  d(kT),  for  a controllable  and 
observable  system  with  Inaccessable  states.  An  observer  Is  designed  for  the 
system  as  in  option  39.  Then  the  control  law  used  is 


u(kT)  * K1  x(kT)  + K2  z(kT) 


(A-22) 


The  composite  closed-loop  system  for  the  disturbance  rejector  is 


■ * 

* 

r « 

x(kT+T) 

ASDHAT +BSDMAT*K1 

,BSDMAT*K2,-BSDMAT*K., 

x(kT) 

z(kT+T) 

TM 

, Ip  . o 

z(kT) 

x(kT+T) 

0 

0 ,ASDMAT+L*CMAT 

x(kT) 

. 

• 

. 

RSDMAT 

0 

0 


d(kT)  (A-23) 


Option  42  Design  Sampled-Data  Tracker 

Option  42  Is  identical  to  option  40  except  that  there  Is  a comnand 
input  instead  of  a disturbance  Input.  Also,  In  this  case  the  addition  of 
discrete- time  Integrators  and  comparators  forces  the  system  to  track  the 
p plecewise-constant  command  inputs  in  the  steady-state  If  the  eigenvalues 
of  the  closed-loop  system  are  assigned  within  the  unit  circle. 

The  control  law  designed  for  the  sampled-data  tracker  Is 

u(kT)  - K1  x(kT)  + Kp  z(kT)  (A-24) 
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and  the  composite  system  formed  using  the  control  law  of  Eq.  A-24  is 


x(kT+T) 

z(kT+T) 


ASDMAT+BSDMAT*^ , BDMAT*K 


TM 


x(kT) 

■ 

0 

+ 

z(kT) 

. 

EMAT 

» ■ 

v(kTJ 


(A-25 ) 


where  EMAT  is  a pxp  diagonal  matrix  with  -T  along  the  diagonal  and  v(kT)  is 
the  piecewise  constant  command  input. 

Option  43  Design  Sampled-Data  Tracker  with  Full-Order  Observer 
Option  43  is  used  when  the  discretized  states,  x(kT),  of  the  system  are 
inaccessable  but  are  observable  In  the  output  vector  y(kT).  In  this 
case,  a full -order  observer  is  designed  for  the  system  as  in  option  38. 
Also,  the  digital  control  law  becomes 


u(kT)  = K1  x(kT)  + K2  z(kT) 


(A-26) 


and  the  composite  closed-loop  system  is 


x(kT+T) 

z(kT+T) 

3 

x(kT+T) 

• i 

ASDMAT+BSDMAT*K1  .BSDMAT^.-BSDMAT*^ 


TM 

0 


,ASDMAT+L*CMAl 


x(kT) 

» ■ 

0 

z(kT) 

+ 

EMAT 

x(kT) 

0 

• ■ 

v(kT) 


Options  44  and  45  Reserve  Option 

2.4  Simulation  Options 

In  order  to  evaluate  the  discrete  controllers  designed  in  options 
28  thru  43,  it  is  necessary  to  perform  a continuous- time 
simulation  on  the  closed-loop  system  with  the  control  input,  u(kT),  held 
constant  over  each  sampling  period  of  the  digital  controller.  Since  it 
is  necessary  to  determine  the  system  state  responses  and  outputs  between 
sampling  times,  options  47  thru  52  discretize  the  continuous  system  with 
a sampling  time  less  than  or  equal  to  T/12,  where  T is  the  sampling  time 
of  the  digital  controller.  Then  the  following  recursive  formula  is  used 

x(kTs  + Ts)  + AS0MATTs  x(kT$)  + BSDMATTsu(kTs)  + RSDMATT$  d(kT$)  (A-28) 

where  T T/12  and  u(kT  ) and  d(kT  ) are  held  constant  over  the  sampling 
Interval  T.  Options  53  to  58  are  discrete  simulations  and  give  the  state 
and  output  responses  only  at  the  sampling  Interval,  T.  Thus,  depending  on 
the  discrete  simulation  choosen,  the  recursive  formula  given  by  one  of  the 
following  equations  Is  used:  Eq.  A-10;  Eq  A-19;  Eq.  A-20;  Eq  A-23; 

Eq.  A-25;  or  Eq.  A-27. 

Option  46  List  Options 

This  option  lists  options  46-60 

Option  47  Continuous-Time  Simulation  of  Regulator  with  Discrete-Time 
Control  Law 

Option  47  performs  a continuous-time  simulation  of  a regulator  with  a 
discrete-time  control  law  given  by 

u(kT)  « K1  x(kT)  (A-29) 


» 
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Option  48  Continuous-Time  Simulation  of  Regulator  with  Discrete-Time 
Control  Law  and  Full -Order  Observer 
This  option  performs  a continuous-time  simulation  of  a regulator  with 
a discrete-time  full-order  observer  and  a digital  control  law  given  by 

u(kT)  = K1  x(kT)  (A-30) 

Option  49  Continuous-Time  Simulation  of  Disturbance  Rejector  with  Discrete- 
Time  Control  Law 

A continuous- time  simulation  of  a disturbance  rejector  with  a discrete- 
time control  law  given  by 

u(kT)  = K1  x(kT)  + ^ z(kT)  (A-31) 

is  performed. 

Option  50  Continuous-Time  Simulation  of  Disturbance  Rejector  with 
Discrete-Time  Control  Law  and  Full-Order  Observer 
This  option  performs  a continuous-time  simulation  of  a disturbance 
rejector  with  a discrete-time  full-order  observer  and  a digital  control 
law  given  by 

u(kT)  = K]  x(kT)  + z(kT)  (A-32) 

Option  51  Continuous-Time  Simulation  of  Tracker  with  Discrete-Time  Controller 
Option  51  performs  a continuous-time  simulation  of  a tracker  with  a 
discrete- time  control  law  given  by 

u(kT)  + K-jX(KT)  + ^ z(kT)  (A-33) 


Option  52  Continuous-Time  Simulation  of  Tracker  with  Discrete-Time 
Controller  and  Full -Order  Observer 

This  option  performs  a continuous-time  simulation  of  a tracker  with  a 
discrete-time  full-order  observer  and  a digital  control  law  given  by 


(A- 34) 


u(kT)  * 1^  x(kT)  + Kg  z(kT) 

Option  53  Discrete-Simulation  of  Regulator 

This  option  uses  the  recursive  formula  of  Eq.  A-10  to  perform  a 
discrete  simulation  of  a regulator  with  a digital  control  law  given 
by  Eq.  A-29. 

Option  54  Discrete  Simulation  of  Regulator  with  Full-Order  Observer 
A discrete  simulation  of  a regulator  with  a full-order  observer  is 
performed,  using  the  recursive  formula  given  by  Eq.  A-17  and  a digital 
control  law  given  by  Eq.  A-30. 

Option  55  Discrete  Simulation  of  Disturbance  Rejector 

Option  55  performs  a discrete  simulation  of  a disturbance  rejector,  using 
the  recursive  formula  given  by  Eq.  A-20  and  a digital  control  law  given  by 
Eq.  A-31. 

Option  56  Discrete  Simulation  of  Disturbance  Rejector  with  Full-Order 
Observer 

This  option  performs  a discrete  simulation  of  a disturbance  rejector 
with  a full -order  observer,  using  the  recursive  formula  given  by  Eq.  A-23 
and  a digital  control  law  described  by  Eq.  A-32. 

Option  57  Discrete  Simulation  of  Tracker 

A discrete  simulation  is  performed  on  a tracker  using  the  recursive 
formula  described  by  Eq.  A-25  and  a digital  control  law  given  by  Eq.  A-33. 

Option  58  Discrete  Simulation  of  Tracker  with  Full-Order  Observer 
This  option  performs  a discrete  simulation  on  a tracker  with  full- 
order  observer,  using  the  recursive  formula  described  by  Eq.  A-27  and  a 
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control  law  described  by  Eq.  A-26. 


| 
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Option  59  Calcomp  Plot  of  Simulation 

This  option  allows  the  user  to  get  Calcomp  plots  of  the  system  state 
or  output  responses  for  any  of  the  above  simulations.  The  Calcomp 
plots  are  stored  in  a local  file  called  PLOT  which  the  user  can  later 
dispose  to  the  plotter  after  CESA  has  been  terminated. 

The  user  is  first  asked  which  simulation  option  is  to  be  plotted. 

Then,  after  that  simulation  has  been  completed  and  the  responses 
printed  out,  the  user  is  asked  whether  he  wants  the  outputs  plotted, 
the  states  plotted,  the  command  inputs  plotted  or  the  disturbance  inputs 
plotted.  After  the  user  has  selected  the  type  of  response  to  be  plotted, 
he  is  asked  the  number  of  the  state,  output,  command  or  disturbance  to 
be  plotted.  At  this  point  the  user  Is  given  the  opportunity  to  control 
the  physical  size  of  the  plot.  Normally  with  FACTORS , the  plot  is 
6x9  inches.  Making  FACTOR  * 2 double  the  size  of  the  plot.  Setting  FACTOR 
* 0.5  reduces  the  plot  to  half  size.  Finally  the  user  is  asked  to 
Input  a plot  title  to  be  drawn  above  the  plot.  The  title  may  be  up  to 
50  characters  long.  A guide  is  printed  to  aid  the  user  as  shown  below. 

!>(-n ENTER  TITLE  <30  CHARACTERS  MAX) )< 

>;£X4MPLE  TRACKER fOU?fUT  #1  RESPONSE 
EXAMPLE  TRACKER  OUTPUT  *1  RESPONSE 

The  program  then  allows  the  user  to  plot  other  responses  or  to  terminate 
the  option.  An  example  plot  with  FACTOR*!  is  shown  in  Fig.  A-2  (see  page  152). 
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2.4  Miscellaneous  Options 

This  option  allows  the  user  to  store  all  the  system  matrices  on  the 
file  MEMORY  during  the  option.  This  is  a precautionary  measure  in  case,  for 
some  unexpected  reason,  the  program  terminates  abnormally. 

Option  100  CESA  Information 

Option  100  gives  a brief  description  of  CESA  and  a group  listing  of 
options. 

3.0  Special  Commands 

At  the  present  time  CESA  has  only  two  special  commands;  STOP;  and 
ANSWER.  These  commands  may  be  typed  only  in  the  option  mode. 

STOP 

This  command  automatically  updates  the  local  file  MEMORY,  and 
prints  out  messages  notifying  the  user  of  any  local  files  that  have  been 
created  during  execution  of  the  program.  Then  CESA  is  terminated. 

ANSWER,  OFF 

When  this  command  is  typed  in,  all  information  is  printed  at 
the  user's  terminal. 

ANSWER,  ON 

This  command  causes  all  program  output,  except  prompt  messages  to 
the  user,  to  be  sent  to  a local  file  called  ANSWER.  When  CESA  is  terminated 
this  file  can  be  routed  to  the  line  printer. 


4.0  CESA  OPTIONS 


O 


COMMAND  - ATTACH,  CESA,  ID=AFIT 
COMMAND  - CESA 

WELCOME  TO  CESA  — VERSION  1.0 
TYPE  100  FOP  HELP 
TYPE  STOP  TO  END  PROGRAM 
OPTION  > 100 

CESA  IS  AN  ITERATIVE  COMPUTER-AIDED  DESIGN  PROGRAM  FOR  DIGITAL  & 
CONTINUOUS  CONTROL  LAW  DEVELOPMENT  USING  COMPLETE  EIGENSTRUCTURE 
ASSIGNMENT.  IT  CONTAINS  60  OPTIONS  THAT  ARE  DIVIDED  INTO  4 GROUPS 
AS  FOLLOWS.  OPTIONS  0,  16,  30,  and  46  GIVE  LISTINGS  OF  THE  OPTIONS  IN 
THESE  RESPECTIVE  GROUPS. 

0-15:  MATRIX  INPUT  OPTIONS 

16-29:  STATE-SPACE  ANALYSIS  OPTIONS 

30-45:  CONTROL  LAW  DESIGN  OPTIONS 

46-60:  SIMULATION  OPTIONS 

OPTION>  0,  16,  30,  46 

CO-15)  MATRIX  INPUT  OPTIONS 

0 *LIST  OPTIONS 

1 ^RECOVER  ALL  DATA  FROM  FILE  MEMORY 

2 *AMAT  CONTINUOUS  PLANT  MATRIX 

3 *BMAT  CONTINUOUS  CONTROL  INPUT  MATRIX 

4 *CMAT  SYSTEM  OUTPUT  MATRIX 

5 *RMAT  CONTINUOUS  DISTURBANCE  INPUT  MATRIX 

6 *CDMAT  SYSTEM  INTEGRATED  OUTPUT  MATRIX 

7 *HELP  USER  SET  UP  CONTINUOUS  STATE-SPACE  MODEL  OF  SYSTEM 

8 *ASDMAT  SAMPLED-DATA  SYSTEM  MATRIX 

9 *BSDMAT  SAMPLED-DATA  CONTROL  INPUT  MATRIX 
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10 

♦RSDMAT  SAMPLED-DATA  DISTURBANCE  INPUT  MATRIX 

w 

ii 

♦HELP  USER  SET  UP  SAMPLED-DATA  STATE-SPACE  MODEL  OF  SYSTEM 

12 

♦RESERVED  OPTION 

13 

♦RESERVED  OPTION 

14 

♦RESERVED  OPTION 

15 

♦RESERVED  OPTION 

(16-29) 

STATE-SPACE  ANALYSIS  OPTIONS 

16 

♦LIST  OPTIONS 

17 

♦COMPUTE  EIGENSTRUCTURE  OF  AMAT 

18 

♦CHECK  CONTROLLABILITY  OF  CONTINUOUS  SYSTEM 

19 

♦CHECK  OBSERVABILITY  OF  CONTINUOUS  SYSTEM 

20 

♦COMPUTE  SAMPLED-DATA  SYSTEM 

21 

♦COMPUTE  EIGENSTRUCTURE  OF  ASDMAT 

o 

22 

♦CHECK  CONTROLLABILITY  OF  SAMPLED-DATA  SYSTEM 

23 

♦CHECK  OBSERVABILITY  OF  SAMPLED-DATA  SYSTEM 

24 

♦COMPUTE  CONTROL  INDICES  OF  SAMPLED-DATA  SYSTEM 

25 

♦COMPUTE  OBSERVER  INDICES  OF  SAMPLED-DATA  SYSTEM 

26 

27 

28 

29 

(30-45) 

30 

31 

32 

33 

34 

35 


♦COMPUTE  CONTROLLABILITY  OF  AUGMENTED  SYSTEM 
♦COMPUTE  OBSERVABILITY  OF  AUGMENTED  SYSTEM 
♦COMPUTE  CONTROL  INDICES  OF  AUGMENTED  SYSTEM 
♦COMPUTE  OBSERVER  INDICES  OF  AUGMENTED  SYSTEM 
CONTROL  LAW  DESIGN  OPTIONS 
♦LIST  OPTIONS 

♦DESIGN  CONTINUOUS  REGULATOR 

♦DESIGN  CONTINUOUS  REGULATOR  WITH  FULL-ORDER  OBSERVER 
♦DESIGN  CONTINUOUS  DISTURBANCE  REJECTOR 

♦DESIGN  CONTINUOUS  DISTURBANCE  REJECTOR  WITH  FULL-ORDER  OBSERVER 
♦DESIGN  CONTINUOUS  TRACKER 
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♦DESIGN  CONTINUOUS  TRACKER  WITH  FULL-ORDER  OBSERVER 

♦RESERVED  OPTION 

♦DESIGN  SAMPLED-DATA  REGULATOR 

♦DESIGN  SAMPLED-DATA  REGULATOR  WITH  FULL-ORDER  OBSERVER 
♦DESIGN  SAMPLED-DATA  DISTURBANCE  REJECTOR 

♦DESIGN  SAMPLED-DATA  DISTURBANCE  REJECTOR  WITH  FULL-ORDER  OBSERVER 
♦DESIGN  SAMPLED-DATA  TRACKER 

♦DESIGN  SMAPLED-DATA  TRACKER  WITH  FULL-ORDER  OBSERVER 

♦RESERVED  OPTION 

♦RESERVED  OPTION 

SIMULATION  ROUTINES 

♦LIST  OPTIONS 

♦CONTINUOUS-TIME  SIMULATION  OF  REGULATOR  WITH  DISCRETE-TIME  CONTROL  LAW 

♦CONTINUOUS-TIME  SIMULATION  OF  REGULATOR  WITH  FULL-ORDER  OBSERVER 
AND  DISCRETE-TIME  CONTROL  LAW 

♦CONTINUOUS-TIME  SIMULATION  OF  DISTURBANCE  REJECTOR  WITH  DISCRETE- 
TIME CONTROL  LAW 

♦CONTINUOUS-TIME  SIMULATION  OF  DISTURBANCE  REJECTOR  WITH  DISCRETE-TIME 
CONTROL  LAW  AND  FULL-ORDER  OBSERVER 
♦CONTINUOUS-TIME  SIMULATION  OF  TRACKER  WITH  DISCRETE-TIME  CONTROL  LAW 

♦CONTINUOUS-TIME  SIMULATION  OF  TRACKER  WITH  FULL-ORDER  OBSERVER 
AND  DISCRETE-TIME  CONTROL  LAW 
♦DISCRETE-TIME  SIMULATION  OF  REGULATOR 

♦DISCRETE-TIME  SIMULATION  OF  REGULATOR  WITH  FULL-ORDER  OBSERVER 

♦DISCRETE-TIME  SIMULATION  OF  DISTURBANCE  REJECTOR 

♦DISCRETE-TIME  SIMULATION  OF  DISTURBANCE  REJECTOR  WITH  FULL- 
ORDER  OBSERVER 

♦DISCRETE-TIME  SIMULATION  OF  TRACKER 

♦DISCRETE-TIME  SIMULATION  OF  TRACKER  WITH  FULL-ORDER  OBSERVER 


♦CALCOMP  PLOT  OF  SIMULATION 
♦UPDATE 
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5.0  Tracker  Design 

In  this  section  a digital  control  law  is  developed  for  the  continuous 
system  that  is  entered  in  the  example  for  option  7.  Options  20,  42, 

59,  and  57  are  used  to  design  the  tracker.  The  following  is  the  actual 
computer  printout  of  the  design  and  simulation.  Also,  the  two  Calcomp 
plots  that  are  set  up  during  the  program  are  included  at  the  end  of  the 
simulation. 
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Appendix  B 


1.  Introduction 

This  manual  provides  documentation  to  assist  In  future  modifications 
and  additions  to  CESA.  This  manual  provides  a description  of  the  overlay 
structure*  subroutines*  and  mass  storage  files  of  CESA  that  Is  essential 
to  any  further  alteration  of  the  program. 

The  reader  Is  assumed  to  have  a working  knowledge  of  the  external 
operations  of  this  program.  A study  of  the  Users  Manual  for  CESA 
should  provide  the  necessary  background  to  understand  the  Internal  parts 
of  the  program  presented  In  this  manual. 

2.  Description  of  Overall  Structure 

The  purpose  of  creating  CESA  Is  to  provide  the  user  with  an  Iterative 
design  tool  for  the  design  of  control  laws  for  regulators,  disturbance 
rejectors*  and  trackers.  Therefore,  It  Is  necessary  that  the  program 
have  the  capabilities  of  designing  the  control  laws  and  evaluating  the 
system  responses.  Also,  the  user  must  have  the  ability  to  redesign  an 
unacceptable  control  law  without  having  to  completely  reenter  all 
pertlnant  data.  Thus,  It  Is  necessary  to  store  all  Information  and  for 
the  program  to  be  Interactive. 

In  order  toe  CESA  to  fullflll  these  requirements.  It  Is  necessary  to 
use  special  space  saving  techniques  so  that  CESA  can  fit  within  the  60kg 
of  memory  that  Is  available  on  INTERCOM.  The  two  space  saving  techniques 
employed  In  CESA's  development  are:  overlays  and  mass  storage. 

2.1  How  Overlays  Work 

Overlays  are  used  to  reduce  the  storage  required  for  a large  program 
by  splitting  the  program  Into  modules.  Each  module  Is  a separate 
program  and  all  the  modules  are  linked  together  through  a main  executive 


overlay.  The  main  executive  overlay  Is  always  In  core  and  controls  the 
program  segments  called  the  primary  overlays.  When  a primary  overlay 
is  needed,  it  is  called  into  central  memory  by  the  main  executive  overlay. 
Only  one  primary  overlay  may  be  in  central  memory  with  the  executive. 
However,  each  primary  overlay  may  control  a group  of  overlays  called 
secondary  overlays.  It  is  possible  to  have  the  executive  overlay,  one 
primary,  and  one  secondary  overlay  in  core  at  one  time. 

Data  is  passed  between  overlays  through  common  statements  or  mass 
storage  files.  In  order  for  data  stored  in  a common  to  remain  accessable 
throughout  the  program  execution,  the  common  statement  must  be  in  the 
main  executive  overlay. 

The  COC  FORTRAN  REFERENCE  MANUAL  and  the  CPC  LOADER  REFERENCE  MANUAL 
contain  more  detailed  information  on  overlays. 

2.2  CESA's  Overlay  Structure 

One  executive  overlay,  14  primary  overlays,  and  16  secondary  overlays 
are  used  In  CESA. 

The  main  executive  overlay  initializes  and  stores  the  data  in  common 
blocks.  Whenever  an  overlay  Is  needed.  It  Is  called  by  the  main  executive 
overlay.  A decision  making  routine  In  the  main  executive  overlay  decides 
when  a particular  primary  overlay  Is  to  be  called.  This  overlay  is 
discussed  further  in  Section. 3.. 

The  primary  overlays,  secondary  overlays,  and  associated  subroutines 
perform  all  the  actual  calculations  In  CESA.  Each  primary  overlay  is 
responsible  for  a main  operation  that  may  require  calling  up  secondary 
overlays  and  subroutines. to  complete  the  operation.  All  subroutines  are 
packed  behind  their  associated  overlay  so  that  they  take  up  memory  only 
when  their  associated  overlay  is  in  main  core.  Primary  and  secondary 
overlays  are  discussed  In  Section  4 and  Subroutines  are  discussed  in 
Section  5. 
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2.3  CESA's  Random  Access  Mass  Storage  Files 

CESA  uses  two  random -access  mass  storage  files:  MEMAUX;  and  SIMMAS. 
These  random-access  files  allow  the  program,  through  the  use  of  special 
read  and  write  statements,  to  store  data  anywhere  on  that  file  and  then. 
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at  a latter  time,  to  recover  the  data. 

The  MEMAUX  file  stores  all  matrices  used  In  the  program  other  than 
the  unaugmented  continuous  and  discrete  system  matrices.  The  random- 
access  files  use  an  Index  to  keep  track  of  the  storage  spaces  used  on  the 
file.  Table  B1  lists  the  name  and  the  index  number  for  the  matrices 
stored  on  the  MEMAUX  file.  Also  listed  in  Table  B1  is  the  location  on 
the  one-dimensional  index  array,  INMASS(I),  that  keeps  track  f the  matrix 
location  and  the  locations  on  the  two-dimensional  array,  MA  T,J), 
that  keeps  track  of  the  current  size  of  the  matrices. 

The  SIMMAS  file  Is  used  to  store  all  simulation  responses.  SIMMIN(I) 
Is  the  Index  array  used  to  keep  track  of  the  location  of  the  different 
system  responses.  Table  B2  gives  the  location  spaces  for  the  system 
responses. 

TABLE  B2 

System  Response  Location  On  SIMMAS  File 


INDEX  NUMBER 

MATRIX 

RESPONSE 

1 

Y 

States 

SIMMINI 

rrr 

2 

YC 

Outputs 

SIMMINI 

2) 

3 

YCDT 

Integrated  Outputs 

SIMMINI 

3) 

4 

YD 

Disturbances 

SIMMIN(4) 

5 

YV 

Commands 

SIMMINC5) 

6 

ROE 

Observer  Error 

SIMMIN 

[6) 

7 

ROB 

Observer  States 

SIMMIN 

(7) 

8 

TI 

Integrator 

SIMMIN 

8) 

9 

1 ^ 

R 

Controls 

SIMMIN 

;9)  - 

More  Information  on  random-access  mass  storage  files  may  be  obtained 


from  the  CPC  FORTRAN  REFERENCE  MANUAL. 
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3.0  Description  of  CESA's  Main  Executive  Overlay 

CESA's  main  executive  overlay  controls  the  calling  of  primary  overlays  and 
therefore,  also  controls  the  flow  of  data  to  these  routines.  Also, 

CESA's  common  storage  locations  are  in  this  main  overlay.  Thus,  the  main 
overlay  is  the  most  important  overlay  In  CESA  and  should  be  thoroughly 
understood  before  any  major  modifications  to  CESA  are  performed. 

3.1  The  Program  Statement 

The  program  statement  is  important  because  it  controls  how  much  buffer 
length  is  allocated  to  reading  and  writlnn  from  the  specified  tapes  used 
In  the  program.  Thus,  if  additional  memory  is  needed  somewhere  else  in  the 
program,  the  buffer  Tength  for  different  tapes  can  be  reduced.  However,  there 
Is  a tradeoff  between  reducing  the  buffer  length  and  Increasing  Input/Output 
time  (T/O).  For  example,  the  Input/Output  Buffer  length  for  the  ANSWER  file 
is  set  to  500g  since  there  is  a large  amount  of  Information  sent  to  this 
file,  while  the  seldom  used  DOODLE  file's  Buffer  length  is  set  to  100^. 

Another  point  to  remember  is  that  there  are  no  restrictions  on  the  amount  of 
I/O  time  that  is  used  on  INTERCOM.  Therefore,  the  only  penalty  Incurred 
for  Increasing  I/O  time  Is  a slow  program  response.  Additional  Information 
on  program  statements  is  available  In  the  CPC  FORTRAN  REFERENCE  MANUAL. 

3.2  CESA's  Common  Data  Base 

Labeled  commons  are  used  in  CESA  for  communication  between  overlays. 
However,  any  Information  or  !data  that  must  be  stored  In  a labeled 
common  for  the  duration  of  the  program,  must  be  stored  In  a labeled  common 
In  the  executive  overlay.  Otherwise  it  will  be  lost  when  the  particular 
overlay  it  originates  in  is  returned  from  central  memory. 


3.3  Description  of  the  Main  Executive  Overlay 

CESA's  main  executive  overlay  controls  the  selection  of  primary  overlays. 
This  section  is  intended  to  describe  how  CESA  performs  this  selection  and 
other  miscellaneous  operations. 

Upon  execution  of  CESA  all  variables. located  In  the  common  data  base  are 
initialized.  Also,  the  two  mass  storage  files  MEMAUX  and  MEMSIM  are 
opened.  CESA,  then  enters  the  option  mode  at  statement  11110  and  calls 
program  READER  (overlay  6)  to  receive  and  Interpret  the  user's  input. 

READER  compiles  the  user's  input,  stores  the  commands  in  coded  form  in 

the  array  MCOMM,  sets  the  MPT  pointer  to  the  first  command,  and  returns  control 

to  the  executive  overlay  at  statement  number  11111. 

Statement  11111  is  the  main  control  node  for  the  program.  All  the 
primary  overlays  return  control  to  this  point.  From  this  node,  the  program 
determines  where  control  should  be  transferred  next.  CESA  uses  a three 
level  sorter  to  accomplish  this.  A detailed  description  of  this  three  level 
sorter  may  be  found  in  Ref.  9. 

4.  Description  of  CESA's  Overlays. 

The  source  listing  of  CESA's  primary  overlays  contains  comment  statements 
which  are  Intended  to  make  the  programs  self-explanatory  to  the  user. 

This  section  contains  an  overview  of  each  of  the  overlays. 

4.1  Overlay  (1 ,0)  — Program  Update 

Program  UPDATE  performs  two  main  functions:  writes  Information  contained 
In  CESA's  common  data  base  to  the  local  file  MEMORY,  and  reads  Information 
on  the  MEMORY  file  into  the  common  data  base. 

Upon  selection  of  option  60,  UPDATE  is  called  by  the  main  executive 
overlay.  The  MEMORY  file  Is  rewound  and  pertalnent  variables  that  are 
stored  in  CESA's  common  data  base  are  written  to  the  file.  A number  code, 

979,  is  written  to  the  tape  so  that  the  MEMORY  file  can  be  identified. 
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Control  Is  returned  to  the  main  executive  overlay  when  this  task  is  completed. 

A reverse  procedure  takes  place  when  option  1 is  selected.  The  MEMORY 
file  is  rewound  and  tested  for  the  proper  number  code,  979.  If  this 
code  is  not  correct,  an  error  message  is  printed,  otherwise  the  information 
on  the  file  is  read  into  the  common  data  base. 

When  the  command  STOP  is  encountered,  UPDATE  performs  option  91  and 
terminates  CESA. 

4.2  Overlay  (2,0)  — Program  MATRIX 


Overlay  (2,0)  is  a short  executive  overlay  which  calls  two  secondary 
overlays  which  are  needed  to  aid  the  user  in  entering  the  system  matrices. 

If  a continuous  system  matrix  is  to  be  entered,  program  MATRIX  calls 
program  MATCON  (Overlay  (2,1)).  Also,  if  a discrete  system  matrix  is  to 
be  entered,  program  MATSD  (Overlay  (2,2))  is  called.  Sections  4.3  and  4.4 
describe  these  two  secondary  overlays. 

4.3  Overlay  (2,1)  — Proqram  MATCON 


Program  MATCON  is  responsible  for  executing  options  that  deal  with  inputting 
continuous  system  matrices.  It  uses  subroutine  MATIN  which  is  discussed 
in  Section  5. 

4.4  Overlay  (2,2)  — Proqram  MATSD 


This  program  is  used  when  an  option  Is  selected  to  input  a sampled- 
data  system  matrix.  MATSD  uses  subroutine  MATIN,  which  is  duscussed  in 


section  5. 

4.5  Overla 


— Program  INDICIE 


Program  INDICIE  is  a short  executive  routine  that  determines  the 
controllability  and  observability  of  a system  and  transfers  control  to  two 
secondary  overlays  Is  the  control  or  observer  Indices  of  a system  are  to 
be  calculated.  Program  TRANSFAB  (Overlay  (4,1))  is  called  to  determine  the 


control  indices  of  a system,  and  program  TRANSFAC  (Overlay  (4,2))  Is  called 
to  determine  the  observer  indices  of  a system. 

4.6  Overlay  (4,1)  — Program  TRANFAB 
Overlay  (4,1)  calculates  the  control  indices  of  a system  by  transforming 

the  system  Into  a canonical  form  and  then  counting  the  size  of  the  system's 
sub-blocks.  Also,  the  determinant  of  the  controllability  matrix  is  determined. 

Program  TRANFAB  uses  subroutines  copy,  MATECHO,  EIGRF,  SETQUP,  UREFORM, 

APLFORM,  BRUNFORM,  SETQDWN,  IDENT.HHPY,  AND  LINV3F  which  are  all  discussed 
In  Section  5. 

4.7  Overlay  (4,2)  — Program  TRANFAC 
This  program  evaluates  the  observer  Indices  of  a system.  TRANFAC  uses 

subroutines  COPY,  TRANPOS.  SETQUP,  UREFORM,  APLFORM,  BRUNFRM,  SETQDWN,  AND 
MATECHO. 

4.8  Overlay  (5,0)  — Program  SIMU 
Program  SIMU  Is  an  executive  routine  that  sets  up  both  the  discrete 

and  continuous  -time  simulations  for  the  designed  closed*loop  system  with 
digital  control  laws.  It  calls  DISSIM  (Overlay  (5,1))  to  perform  the  discrete 
simulation  and  it  calls  CONSIM  (Overlay  (5,2))  to  perform  the  continuous-  , 

time  simulation.  After  the  simulation  has  been  completed  and  control 

1 

returned  to  SIMU,  program  PRNTSIM  Is  called  to  list  the 
results. 

j 

4.9  Overlay  (5,1)  — Program  DISSIM 

Overlay  (5,1)  performs  the  discrete  simulation  using  a sampled-data 
recursive  formula  representation  of  the  closed-loop  system. 


° 

I 


4.10  Overlay  (5,21  — Program  CONSIM 

This  program  performs  a continuous-time  simulation  on  the  closed-loop 
system  as  described  in  the  USER'S  MANUAL  for  CESA. 

4.11  Overlay  (.5,3)  — Program  PRNTSIM 

PRNTSIM  lists  all  the  results  of  the  closed-loop  simulation. 

4.12  Overlay  (.6,0)  --  Program  READER 

Program  READER  verifies  and  codes  all  input  that  Is  entered  in  the 
option  mode.  Reference  9 gives  a detailed  description  of  READER. 


4.13  Overlay  (7.0] 


Program  DECODER 


Program  DECODER  is  used  to  decode  the  MCOMM  and  DATM  arrays  set  up  by 
READER.  Ref.  9 gives  a detailed  description  of  this  routine. 


This  program  is  not  being  used  presently  in  CESA. 


4.15  Overlay  (9,0) 


Program  LISTOPT 


Program  LISTOPT  contains  the  print  statements  that  list  the  options  of 
CESA. 

I 

4.16  Overlay  ClOiO)  — Program  SDREG 

Program  SDREG  is  the  main  routine  used  to  perform  the  sampled-data 
regulator  design  options.  It  calls  on  the  two  secondary  overlays. 

NSAB1  and  DISCLS1 , to  calculate  the  control  law  and  to  set  up  the  closed-loop 
system. 

4.17  Overlay  00,1)  — Program  NSAB1 

This  program  contains  the  coded  algorithm  for  entire  eigenstructure 
assignment.  It  allows  the  user  flexibility  In  assigning  eigenvalues  and 
the  associated  eigenvectors.  This  program  is  called  to  calculate  the  state 
feedback  gain  matrix  and  also  to  calculate  the  observer  gain  matrix. 
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4.18  Overlay  Cl 0.2)  — Program  DISCLS1 

Overlay  Cl 0 , 2 ) sets  up  the  closed-loop  system  for  the  designed  control 

law. 

4.19  Overlay  (11,0)  — Program  SDDREJ 

Program  SDDREJ  is  the  executive  routine  used  to  calculate  the  digital 
control  law  for  the  disturbance-rejector  system.  It  calls  program  NSAB  and 
program  DISCLS  to  determine  the  digital  control  law  and  set  up  the  closed-loop 
system. 

4.20  Overlay  (.H  .1)  — Program  NSAB 

Program  H5AB  4s  identical  to  program  NSAB1. 

4.21  Overlay  (11,2)  — Program  DISCLS 

Overlay  (11,2)  Is  identical  to  program  DISCLS1. 

4.22  Overlay  (12,0)  — Program  SDTRAC 

Program  SDTRAC  is  the  executive  routine  used  to  calculate  the  digital 
control  law  for  the  tracker  system.  Program  NSAB2  and  program  DISCLS  2 are 
called  by  SDTRAC  to  determine  the  digital  control  law  and  to  form  the  closed- 
loop  system. 

4.23  Overlay  (12,1)  — Program  NSAB2 . 

This  program  Is  identical  to  program  NSAB1. 

4.24  Overlay  (12.2)  — Program  DISCLS2 

Overlay  (12,2)  Is  identical  to  program  DISCLS1. 

4.25  Overlay  (13,0)  - Program  DISCRETE 

Overlay  (13,0)  is  the  overlay  that  checks  which  type  of  digital  control 
law  is  to  be  calculated  and  calls  the  appropriate  overlay  accordingly. 

Also,  this  routine  checks  whether  the  sampled-data  system  has  been  entered 
and,  if  it  has  not  been  entered,  it  calls  program  SAMPLE. 


4.26  Overlay  Cl 5,0)  --  Program  SAMPLE 

Program  SAMPLE  is  the  primary  overlay  that  is  called  to  discretize  the 
continuous  system.  It  calls  the  two  secondary  overlays,  EIGVAY  and  EIGSTR, 
to  form  the  sampled-data  system.  Also,  this  overlay  performs  the 
options  17  and  21. 

4.27  Overlay  (15,1)  --  Program  EIGVAV 


Overlay  05,0)  calculates  the  eigenvalues  and  eigenvectors  of  the  plant 


matrix. 


The  secondary  overlay  CALPLOT  uses  Cal  comp  routines  to  plot  the 
data  sent  to  it  by  program  PLOTTS. 

5.0  Description  of  Subprograms 


Overlay  (15,2)  sits  up  the  inverse  modal  matrix. 

4.29  Overlay  (16,0)  — Program  PLOTTS 

Program  PLOTTS  Is  used  when  option  59  is  selected.  It  sets  up  the 
arrays  to  be  plotted  from  the  simulation  responses.  Once  it  sets  up 
the  arrays  It  calls  the  secondary  overlay,  CALPLOT,  to  plot  the 
responses. 

) — Program  CALPLOT 


An  inJdepth  description  of  many  of  the  subroutines  used  In  CESA  may  be 
found  in  Refs.  3,  9 and  15. TABLE  B-l  lists  these  subroutines  and  the 
references  where  they  may  be  found.  Descriptions  of  the  subroutines  used 
in  CESA  that  are  not  found  In  the  preceding  references  are  found  In  the 
following  sections. 


TABLE  B-l 


Subroutine  Reference  List 


SUBROUTINE 

REFERENCE 

SUBROUTINE 

REFERENCE 

MADD 

TOTAL CRef  9) 

SETQDWN 

' FORTRAC(Ref  3) 

INVERT 

TOTAL.(R€f"9) 

CINDS 

FORTRAC(Ref  3) 

READN 

TOTAL (Ref  9) 

APLFORM 

FORTRAC(  Ref'S) 

TRANPOS 

TOTAL (Ref  9) 

UREFORM 

FORTRAC(Ref  3) 

MMPY 

TOTAL (Ref  9) 

SETQUP 

FORTRAC(Ref  3) 

MATECHO 

TOTAL (Ref  9) 

FULLOB 

FORTRAC(Ref  3) 

IDENT 

TOTAL (Ref  9) 

AUGMAT 

FORTRAC(Ref  3) 

MATIN 

TOTAL (Ref  9) 

M3CMAT 

FORTRAC(Ref  3) 

LISTER 

TOTAL (Ref  9) 

ADDMATX 

FORTRAC(Ref  3) 

TVIODV 

TOTAL (Ref  9) 

DE01CKF 

FORTRAC ( Ref  3) 

MODIFY 

TOTSL ( Ref  9) 

SUBMATX 

FORTRAC(Ref  3) 

TITLES 

TOTAL (Ref  9) 

STIMEST 

FORTRAC (Ref  3) 

CONVERT 

TOTAL (Ref  9) 

STSTATE 

FORTRAC (Ref  3) 

SUB1 

FORTRAC(Ref  3) 

SAMRES 

FORTRAC (Ref  3) 

SUB2 

FORTRAC(Ref  3) 

EIGRF 

IMSL(Ref  18) 

UFORM 

FORTRAC(Ref  3) 

LINV3F 

IMSL(Ref  18) 

CHEROW 

FORTRAC(Ref  3) 

BALBNC 

EISPACK(Ref  15) 

ACHEROW 

FORTRAC (Ref  3) 

ELMHES 

EISPACK(Ref  15) 

SUBIOOO 

FORTRAC(Ref  3) 

ELTRAN 

EISPACK(Ref  15) 

SUB2000 

FORTRAC(Ref  3) 

HQR2 

EISPACK(Ref  15) 

SUB3000 

FQRTRAC(Ref  3) 

BALBANC 

EISPACK(Ref  15) 

BRUNFRM 

FORTRAC(Ref  3) 

5.1  Subroutine  AVGAIN 


This  subroutine  determines  the  average  absolute  magnitude  of  the 
feedback  gain  matrix 

Calling  Sequence 

CALL  AVGAIN  (XK,M2,M1 .NKX.MKX) 

Definition  of  Symbols 
XK  3 Feedback  gain  matrix 
M2  3 Number  of  rows  of  XK. 

Ml  3 Number  of  columns  of  XK 

NKX  3 Number  of  rows  of  XK  dimensioned  in  main  routine 
MKX  3 Number  of  columns  of  XK  dimensioned  in  main  routine 

5.2  Subroutine  FACTORL 

This  subroutine  calculates  the f&cfcOrtal  of  a number. 

Calling  Sequence 

CALL  FACTORL  (N,  NFACT , IFAIL) 

Definition  of  Symbols 

N 3 Integer  for  which  factorial  is  to  be  calculated. 

NFACT  3 N1 


IFAIL  3 Error  Flag  . IFAIL  equals  1 if  N is  within 
proper  ranges.  IFAIL  equals  2 if  N ^ 0 

5.3  Subroutine  EXPAT 


and  the 
system. 


The  subroutine  EXPAT  calculates  the  state  transition  matrix  e 
At 


AT 


j: 


e dt.  This  routine  is  used  in  discretizing  a continuous  li 


Calling  Sequence 

CALL  EXPAT  (N,A,DEL,EA,EAINT, IFAIL) 


Appendix  C 

Matrices  for  the  Three  Different  Flight  Conditions 

This  Appendix  contains  some  of  the  continuous  discretized  system 
matrices  for  the  regulator  and  tracker  design  of  Chapter  V according 
to  the  following  difference  equation,  for  a sampling  time  T = 2 seconds 


x(2k+2)  = F x(2kl  + G u(2k) 


(C-l) 


For  both  the  low  altitude  cruise  condition  regulator  and  tracker  the 
continuous  plant  from  Eq.  55  is 


1.6754,-  .755E-2,-.250E-3,  0 3.954  , 1 .033  , 0 

544.9  ,-1.135  rO.1179  , 0 ,-30.77  ,13.5212,  0 

4.374  , .02961  .-.917E-2.-32.18,-  .2470,-2.592  ,-2.6 


1 , 0 


[ 0 ,0  ,0  ,0,0 
and  the  discretized  plant  matrices  are 


,0  , 0 ,0  ,0 

, 0 ,-10  ,0  ,0 

,0,0  ,-4  ,0 

,0,0  ,0  -2 


(C-2) 


- .0325,  .0002,  .0001,-  .0035,  .0189  ,-  .0180,-  .0001 

-11.98  ,-.0477, -.0555,  3.364  ,4.399  ,-1.203  , .0749 

-11.23  , .0632,  .9777,-63.68  ,3.983  ,-2.928  ,-1.255 

.1789, -.0013,  .0001,  .9984, -.0685  ,0  ,0 

0 ,0  ,0  ,0  , .206E-08,  0 , 0 


,0  ,0  ,0  ,0 

,0  ,0  ,0  ,0 


, .0003,  0 

, 0 , .0183 


^ . r‘-  ■*  »• 


(C-3) 


Using  the  stability  derivatives  from  Table  II  and  Eq.  55,  the  continuous 
plant  matrix  for  both  the  regulator  and  tracker  for  the  medium  altitude 
cruise  condition  becomes 


The  discretized  plant  matrices  for  both  the  regulator  and  tracker  for 
the  medium  altitude  cruise  condition  are 


0780,-  .0001,0,0,0 


The  augmented  system  for  the  medium  altitude  condition  Is 


- .0944,  .0003,  .0001,-  .0105,  .0430  ,-  .0303,  - .0002,0  fi  ,0 

-34.01  , -.1110,  .0545,-  1.827  ,10.19  ,-2.547  ,-.0493.0  A ,0 

-11.30  , .0738,  .9803,-63.83  , 3.606  ,-2.320  ,-1.016  ,0  . A ,0 

.1328,  -.0014,  .0003,  .9904,-  .0471  ,0  ,0  ,0  A ,0 

0 ,0  ,0  ,0  , .206E-8,  0 ,0  ,0  A ,0 

0 , 0 ,0  , 0 , 0 , .0003,  0 ,0  ,0  ,0 


,1320  , 0 


, 0 , .0183,0  ,0  .0 

,0  ,0  .1.0,0  ,0 

, 0 ,0  ,0  .1.0,0 

, 0 ,2  .0  A ,1.0 


■ - 


[F+GK] 


(C-8) 


- .1326,  .0008 

, .217E-5 

,-  .5224, 

.0563, 

- .0410.-.152E-4 

-68.20  , .3153 

,-.0620 

,-290.8  , 

25.38  , 

-14.20  , .0643 

- 1.014  , .0050 

, .2293 

,-  7.118  , 

.3158, 

- .2971 ,-.2770 

[F+GK]* 

.0091 ,-.0001 

.-.299E-4 

, .0386, 

- .0050, 

.0053,  .367E-4 

- .2686,  . 952E-3 ,- . 2486E-3 

,-  .9896, 

.0999, 

- .0527,  .276E-3 

- 1.540  , .0070 

,-.0014 

,-  6.492  , 

.5739, 

- .3202,  .0014 

- .4223,  .0030 

, .1907 

,-  4.170  , 

.1216, 

-1923  ,-.2302 

and  the  closed-loop  matrix  for  the  medlam  altitude  regulator  Is 

- .1831,. 0009 

.-.925E-4, 

.5445, 

.0712,- 

.0438  , .349E-4 

-112.6  ,.5014 

,-.0698  ,- 

512.7  ,35 

.58  ,-15.79  , .0445 

- .7602,. 0052 

, .2303  ,- 

4.368  ,- 

.2353,- 

.2391  ,-.2244 

[F+GK]* 

- .0319,. 0001 

.-.234E-4,- 

.1856, 

.0076,- 

.956E-3,  .201E-4 

- .2919..639E-3 

.-.148E-3,- 

1.131  , 

.0837,- 

.0259  , .120E-3 

- 1.891  ,.0075 

,-.0015  ,- 

8.908  , 

.5641,- 

.2159  , .0011 

- .5140,. 0049 

, .2368  ,- 

2.912  , 

.1887,- 

.2538  ,-.2308 

In  the  case  of  the  tracker  designs,  the  closed-loop  matrix  for  the  low 
altitude  tracking  system  Is 


I 


(C-9) 
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1607,  .0019,  .266E-4,-  1.169 


» 


.0675,-  .0762.-.183E-4, 


» 


-85.87  , .9554,  .0218  ,-660.1  ,32.43  ,-33.0  , .0481  , 

- .6766,  .0077, -.7164  ,-  5.899  , .2413,-  .0802, -.0637  , 

- .1109,  .0023,  .108E-3,-  1.428  , .0412,-  .0282,  .591E-4, 

[F+GK]  * .0283,  .0016.-.309E-3,-  .8049,-  .0089,  - .2810,  .412E-4, 

- .5901,  .0156, -.0012  ,-  9.305  , .2308,-  1.358  , .473E-3, 


.8975, -.01 57,  .4286  , 4.290  , .2690,  .8667, -.2833  , 

0 r2.0  .0  .1090  . 0 . 0 .0 


Appendix  D 
C*  Criterion 


In  order  to  determine  a meaningful  performance  criterion  for  a particular 
control  system  it  is  necessary  to  determine  the  types  and  levels  of  handling 
qualities  that  are  required  Ey  a pilot  to  accomplish  a mission.  A particular 
criterion  that  was  developed  to  specify  short  period  handling  qualities 
in  terms  of  aircraft  parameters  that  are  familiar  to  the  pilot  is  the  C* 
criterion  (Ref:16). 

A study  was  done  in  this  thesis  to  determine  if  a deadbeat  controller 
could  be  designed  to  track  a C*  commanded  input.  However  the  augmented 
tracking  system  had  a low  degree  of  controllability  and  the  C*  tracker 
design  was  dropped.  Time  limitation  prevented  an  examination  of  the 
multivariable  zeros  of  the  system  to  determine  their  location  with  respect 
to  the  unit  circle.  Tn  order  to  provide  some  background  for  this  problem, 
the  C*  equations  used  for  the  three  sets  of  trackers  are  included.  There 
was  no  speed  control  for  the  C*  tracker  designs. 

The  C*  equation  used  for  all  three  tracker  designs  was 

c*  = *"zp  + k2  q (D"] ) 

where  the  normal  acceleration  of  the  pilot, 
finzp  in  units  of  g's  is 

Anzp  - (w  - UQ0  - lx  V )/32 .18  , (D-2) 

fcx  is  the  distance  of  the  pilot  from  the  aircraft's  center  of  gravity, 
q is  the  aircraft's  pitch  rate,  and 
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where  C-  Is.  a dimensional  constant  given  in  Table  D1  and  V is  the 
i co 

crossover  velocity  as  defined  in  Ref.  16.  Using  Eqs.  D-l  and  D-2  and 
Table  D1  with 


= 350  ft/sec 


the  C*  equation  for  the  three  flight  condition  was:  for  the  landing 
condition 


C*  = -.2809u-.3243w+27.93<ie+22. 72$  +47.456 


for  the  low  altitude  cruise  condition 


C*  * - . 003298u- . 0242W+4 . 81 86e-l . 0896  +2.6630 


and  for  the  medium  altitude  cruise  condition 


C*  » - . 00247u- . 01 534w+4 . 272de- . 83296  +1.93630 
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